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Chapter  1 

Executive  Summary 


This  report  summarizes  the  work  on  “Joint  Space-Time  Coded  Modulation  and  Channel 
Coding  over  Fading  Channels  with  Cochannel  Interference  ”  carried  out  at  the  NJIT  from 
July  1999  to  August  2002.  The  main  contribution  of  this  work  was  to  develop  algorithms  for 
adaptive  arrays  for  wireless  communications  and  analyze  their  performance  in  fading  channels 
with  cochannel  interference.  The  research  covers  several  topics  that  are  summarized  below. 

We  first  studied  the  performance  of  coherent  detection  scheme  for  adaptive  arrays  over 
wireless  channel.  For  adaptive  arrays,  optimum  combining  (OC)  is  an  excellent  processing 
technique  to  combat  multipath  fading  and  suppress  cochannel  interference.  In  this  report 
we  present  new  methods  to  derive  closed-form  expressions  for  the  exact  bit  error  probability 
(BEP)  for  optimum  combining  ([1],  [2]).  Our  method  differs  from  previous  ones  in  that  it 
starts  from  the  decision  metrics  of  OC  instead  of  the  moment  generating  function.  This 
approach  facilitates  obtaining  closed-form  expressions.  The  final  expression  is  for  multiple 
interferers  and  multiple  reception  branches,  with  the  number  of  interferers  less  than  the 
number  of  reception  branches.  We  assume  that  the  modulation  is  binary  phase  shift  keying 
(BPSK),  and  the  channels  are  independent  Rayleigh  fading  channels.  With  these  expressions 
for  BEP,  evaluating  the  performance  of  OC  is  easy  and  accurate. 

OC  is  a  coherent  detection  scheme.  To  implement  it,  the  following  channel  information 
is  required:  the  channel  gain  (including  amplitude  and  phase)  of  the  desired  signal,  the 
covariance  matrix  of  the  interference  plus  noise.  For  communication  systems  where  channel 
phase  is  very  difficult  or  impossible  to  recover,  OC  is  not  practical  and  noncoherent  detection 
is  necessary.  We  propose  a  noncoherent  detection  scheme  for  adaptive  arrays  ([3],  [4]).  The 
scheme  is  called  multiple  symbol  differential  detection  (MSDD),  in  which  the  channel  gain  of 
the  desired  signal  is  unknown  to  the  receiver,  but  the  covariance  matrix  of  the  interference 
plus  noise  is  known.  The  maximum  likelihood  decision  statistic  is  derived  for  the  detector 
and  its  performance  is  demonstrated  by  analysis  and  simulation.  To  reduce  the  computation 
complexity  of  the  MSDD  decision  statistic,  we  present  a  sub-optimum  decision  feedback 
algorithm  for  iterative  symbol  detection.  This  sub-optimum  algorithm  achieves  performance 
that  is  very  close  to  that  of  optimum  algorithm.  A  closed-form  approximation  to  the  union 
bound  of  pairwise  error  probability  is  shown  to  provide  a  good  approximation  to  the  bit 
error  probability.  We  show  analytically  and  numerically  that  with  an  increasing  observation 
interval,  the  performance  of  this  kind  of  MSDD  approaches  that  of  OC  with  differential 
encoding.  We  also  develop  MSDD  for  another  kind  of  noncoherent  detection  where  the  only 
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required  channel  information  is  the  channel  amplitude  of  the  interference.  It  is  shown  that 
when  the  interference  level  is  high,  this  MSDD  technique  can  achieve  good  performance. 
This  work  was  presented  in  [5]. 
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Chapter  2 


Introduction 


2.1  Motivation 

In  modern  commercial  wireless  communication  systems  such  as  code  division  multiple  access 
(CDMA)  systems  and  time  division  multiple  access  (TDMA)  systems,  the  cellular  concept  is 
widely  applied  to  increase  system  capacity  [6].  Under  this  concept,  the  entire  service  area  is 
divided  into  small  areas  called  cells.  Several  cells  comprise  a  cell  cluster.  A  cell  cluster  uses 
all  the  available  resource.  For  TDMA  the  resource  is  frequency  channels,  while  for  CDMA  the 
resource  is  codes.  For  TDMA  systems,  each  cell  can  use  a  portion  of  the  available  frequency 
channels,  but  cells  in  the  neighborhood  use  different  frequency  channels.  The  cells  that  use 
the  same  frequency  channels  are  at  least  a  cell  away  from  each  other .  In  this  way  the  limited 
precious  frequency  resource  could  be  reused,  hence  the  system  capacity  can  be  increased 
infinitely  (at  least  in  theory),  while  at  the  same  time  the  interference  is  kept  to  a  minimum. 
What  should  be  kept  in  mind  is,  though  the  interference  is  minimized,  it  still  exists  due  to 
the  same  frequency  channel  used  by  different  cells.  Actually  the  interference  becomes  one  of 
the  factors  that  limit  the  performance  of  the  wireless  communication  systems. 

The  cellular  concept  is  illuminated  in  Fig.  2.1.  A*  to  Gj  (i  =0,  1,  ...,  6)  are  cells  that 
form  a  cell  cluster.  Cells  A<  for  i  =0,  1,  ...,  6  use  the  same  frequency  channels,  so  do  cells  B* 
to  Gj.  The  user  sIti  in  cells  A;  (i  =1,  2,  ...,  6)  can  interfere  the  user  s  in  cell  A0. 

Another  factor  that  constraints  the  performance  is  multipath  fading.  In  the  wireless 
environment,  due  to  reflection,  diffraction  and  scattering,  the  transmitted  signal  may  fade 
greatly  and  reach  the  receive  antenna  through  more  than  one  path  (shown  in  Fig.  2.2). 
Since  the  locations  of  the  transmitter,  the  obstacles  and  the  reflectors  are  random,  the 
transmit  paths  are  random  as  well.  The  received  signal  from  a  different  path  may  add  up 
constructively,  or  destructively.  The  total  effect  of  this  summation  is  a  random  attenuation 
of  the  transmitted  signal.  When  the  attenuation  is  deep,  the  received  signal  is  so  weak 
that  the  receiver  won’t  be  able  to  recover  the  transmitted  signal.  To  resolve  this  problem, 
diversity  is  introduced.  By  using  diversity  technique,  several  replicas  of  the  same  information 
signal  are  transmitted  over  independently  fading  channels.  The  probability  that  all  the  signal 
components  reaching  the  receiver  will  fade  simultaneously  is  reduced  considerably. 

Three  of  the  techniques  that  can  achieve  more  than  one  independently  fading  version  of 
the  same  information-bearing  signal  are  [7,  Chapter  14]: 
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Figure  2.1:  The  cellular  concept.  Since  cells  Aj  for  i  =  1,  ...,  6  use  the  same  frequency 
channels  as  cell  A0,  the  user  sjti  (i  =1,  2,  ...,  6)  may  interfere  the  user  s  in  cell  A0 


Figure  2.2:  Multipath.  The  transmitted  signal  reaches  the  receiver  through  path  1  and  path 
2. 
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•  Temporal  diversity:  the  same  information  bearing  signal  is  transmitted  in  more  than 
one  time  slot,  where  the  separation  between  successive  time  slots  equals  or  exceeds  the 
coherence  time  of  the  channel. 

•  Frequency  diversity:  the  same  information  bearing  signal  is  transmitted  on  more 
than  one  carrier  frequency,  where  the  separation  between  successive  carrier  frequencies 
equals  or  exceeds  the  coherence  bandwidth  of  the  channel. 

•  Spatial  diversity:  more  than  one  transmit  and/or  receive  antenna  are  employed.  The 
antennas  are  spaced  sufficiently  far  apart  that  the  multipath  components  in  the  signal 
have  independent  fading. 

Since  spatial  diversity  does  not  require  expansion  of  the  bandwidth,  it  is  desirable  for  a 
bandwidth-limited  system  when  cost  and  size  permitted.  And  as  pointed  out  in  [8],  spatial 
diversity  could  be  used  to  cancel  interference  as  well  as  to  combat  fading.  Capacity  of  a 
system  with  spatial  diversity  has  been  proven  to  increase  with  the  number  of  antennas  [9]. 

It  is  for  these  advantages  that  communication  systems  with  spatial  diversity  has  been  an 
appealing  research  area.  Reception  adversity  has  been  implemented  in  base  stations.  For 
example,  in  the  second  generation  IS-136  TDMA  base  station  [10],  two  receive  antennas 
are  deployed.  Technology  has  been  developed  for  deploying  4  receive  antennas  at  the  base 
station. 

Currently  most  of  the  research  on  spatial  diversity  is  focused  on  space- time  codes  ([11], 
[12],  [13]),  which  employ  transmit  diversity.  While  space-time  codes  can  provide  some  coding 
gain  as  well  as  spatial  diversity,  and  could  be  the  future  application,  this  work  focuses  on 
a  more  practical  problem  for  now:  performance  analysis  of  communication  systems  with 
reception  diversity  but  without  transmit  diversity  and  coding.  To  more  easily  understand 
the  topic,  the  following  presents  the  basic  system  model  used  in  this  work. 


2.2  System  Model 

We  consider  communication  system  with  reception  diversity  but  without  transmit  diversity. 
All  the  signals  are  baseband  signals.  As  shown  in  Fig.  2.3,  there  is  one  transmit  antenna, 
L  reception  branches  and  Nj  interferers  in  the  system.  The  sampled  output  of  the  matched 
filter  for  the  Ath  branch  at  time  k  is  expressed  as 

Nj 

rk,t  =  \/PaCiSk  +  '^2  +  nM»  ,L  (2.1) 

i=l 

where  the  parameters  in  (2.1)  are  defined  as: 

Ps :  power  of  the  desired  signal. 

ci :  channel  gain  of  the  £-th  branch  for  the  desired  signal. 

Sk  :  desired  transmitted  symbol. 

Pj  :  power  of  the  interferers. 

d  i :  channel  gain  of  the  £-th  branch  for  the  z-th  interferer. 

Si,*;  :  interference  signal. 
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Figure  2.3:  Diagram  for  systems  with  L  reception  diversity  branches.  sk  is  the  desired  signal. 
There  could  be  more  than  one  interferer  si>k  (only  one  is  shown  in  the  figure). 


nktt :  complex  white  Gaussian  noise. 

S).  and  Si,k  could  be  multiple  phase  shift  keying  (M-PSK)  symbols,  differential  multiple 
phase  shift  keying  (M-DPSK)  symbols,  or  Gaussian  distributed  signals.  We  will  define  them 
more  specifically  in  later  chapters. 

The  signal  model  in  vector  notation  is 


A h 

y/~PsCSk  "f*  (2.2) 

t=i 

where  r*  =  [r^i,  rkl2,  •  •  •  ,  Tk,L,]T ,  and  the  superscript  T  denotes  vector  transposition;  c,  c{ 
and  n*  are  vectors  that  are  similarly  defined  as  r*,. 

The  channel  gains  c/  and  c*,/  are  assumed  to  be  independently  and  identically  distributed 
(i.i.d.),  zero-mean,  circularly  symmetric,  complex  Gaussian  random  variables  (Rayleigh  fad¬ 
ing))  with  variance  Clt/2  (for  c/)  and  1/2  (for  c^/)  per  dimension.  For  future  use,  we  also 
define  the  fading  matrix  for  the  interferers  as  C/  =  [ci,  C2,  •  •  •  ,  c^;] . 

Define  the  interference  plus  noise  vector  as 


N, 

z*  =  \fPj  ^2  Ci si)k  +  nfc  (2.3) 

»=i 

then  (2.2)  becomes 


rk  =  VPsCSk  +  Zjfe  (2.4) 

In  the  analysis,  we  often  treat  c j  as  constant  vector  first,  then  treat  it  as  random  vector. 
When  we  treat  c /  as  constant  vector,  the  covariance  matrix  of  z*.  is 


Ni 

R z  =  Pi^2  c'ci  +  diag 

i=l 


( J1  2  _2  \ 

V^l  5  a2  j  *  #  ‘  iaL) 


(25) 
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Figure  2.4:  Diagram  of  the  optimum  combining  detector. 


where  the  superscript  H  denotes  the  Hermitian  transposition  and  (of ,  of, • • •  ,  of)  is  the 
power  profile  of  the  noise. 

2.3  Background 

For  wireless  communication  systems  with  reception  diversity,  Optimum  combining  (OC)  is 
a  well-known  approach  to  combat  fading  and  suppress  cochannel  interference.  Its  structure 
is  shown  in  Fig.  2.4. 

In  Fig.  2.4,  the  received  signal  r*  =  [r^i,  r*>2,  •  •  •  ,  rfcj£,]r  is  weighted  and  combined.  The 
weight  vector  w  =  [w1,tn2,  •  •  •  ,  wl]T  is  chosen  in  such  that  the  output  of  combiner  (which 
is  the  input  to  the  detector)  achieves  the  maximum  signal-to-interference-plus-noise  ratio 
(SINR).  The  decision  rule  for  OC  is 

sk  =  arg  max  p  (r*|sfc,  c,  Rz)  (2.6) 

where  p(r*|sfc,  c,Rz)  is  the  probability  of  r*  conditioned  on  sk,  c,  and  R*.  A  simplified 
version  of  this  decision  rule  will  be  shown  in  Chapter  3. 

One  of  our  efforts  is  to  derive  closed-form  expressions  for  bit  error  probability  (BEP) 
for  OC.  Those  kinds  of  expressions  have  only  been  obtained  for  some  special  cases.  Some 
related  work  about  OC  is  summarized  in  Chapter  3. 

OC  is  a  coherent  detection  scheme.  To  construct  the  weight  vector  w,  the  following 
channel  information  is  required:  c,  which  is  the  channel  gain  (include  amplitude  and  phase) 
of  the  desired  signal  and  R2,  the  covariance  matrix  of  the  interference  plus  noise.  For 
communication  systems  where  channel  phases  are  very  difficult  or  impossible  to  recover,  OC 
is  not  practical.  Under  this  circumstance,  noncoherent  detection  scheme  must  be  considered. 

One  kind  of  noncoherent  detection  scheme  is  differential  detection,  in  which  the  transmit¬ 
ted  signals  are  differentially  encoded.  For  conventional  differential  detection,  two  received 
signals  are  used  at  the  observation  interval  to  make  decision  about  the  transmitted  signal. 
The  recovery  of  channel  phase  is  not  required.  The  decision  rule  for  conventional  differential 
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detection  is 


(Slb-i.Slb)  =  arg  max  p(r*_i,rfc|sjfe_i,sfc,Rj , 

3jfe-l,Sfc 


(2.7) 


where  p(rfc_i,rfc|sfc_i,s*;,Rz)  is  the  probability  of  rfc_i,rfc  conditioned  on  sk-i,sk  and  Rz. 
Conventional  differential  detection  suffers  performance  penalty  comparing  to  coherent  de¬ 
tection. 

Multiple  sysmbol  differential  detection  (MSDD)  achieves  performance  between  that  of 
conventional  differential  detection  and  coherent  detection.  In  MSDD,  more  than  two  symbols 
are  used  in  the  observation  interval.  It  was  shown  that  with  the  increase  of  symbols  in  the 
observation  interval,  the  performance  of  differential  detection  can  be  improved  significantly. 
The  decision  rule  for  MSDD  is 


sfc  =  argmaxp  (rjs*,  Rz) ,  (2.8) 

Sfc 

rp 

where  sfc  =  [sk-(K-i),  •  •  •  ,  sk-i,  s*]T  is  a  sequence  of  K  (K  >  2)  symbols,  rfc  =  [rfc_(Ar_i),  •  •  •  ,  rfc] 
is  a  vector  of  all  the  received  signal  in  the  observation  interval,  and  p  (rfc|s*,  RJ  is  the  prob¬ 
ability  of  r*  conditioned  on  s*  and  R2.  Some  related  work  about  MSDD  is  summarized  in 
Chapter  5. 

In  this  work,  we  applied  MSDD  to  communication  systems  with  interference.  The  only 
required  channel  information  for  that  kind  of  MSDD  is  the  covariance  matrix  of  the  interfer¬ 
ence  plus  noise  R*.  By  simulation  and  analysis  results,  we  demonstrate  that  MSDD  could 
achieve  performance  close  to  that  of  OC  with  differential  encoding. 

We  also  develop  MSDD  for  another  kind  of  noncoherent  detection  where  the  only  required 
channel  information  is  the  channel  amplitude  of  the  interference.  Its  decision  rule  is 

(sfc,  s/,fc)  =  arg  max  p  (rjs*,  sI>k,  |c/|) ,  (2.9) 

where  sk  and  r*  are  defined  as  mentioned  before;  s i<k  =  [s/,fc-(K-i),  ■  •  •  ,  S/,*]r  is 

the  sequence  of  interference  symbols;  |c»|  is  the  channel  amplitude  of  the  interference,  and 
p  (ijtlsfc,  S[tk,  |Cj|)  is  the  probability  of  r*  conditioned  on  s*,  s i<k,  and  |Ci|. 

It  is  shown  that  when  the  interference  level  is  high,  this  MSDD  technique  can  achieve 
better  performance  than  detectors  using  optimum  combining  (with  differential  encoding). 


2.4  Outline  of  the  Report 

The  main  topics  of  this  work  are: 

•  Performance  analysis  of  OC. 

•  Derivation  of  the  decision  statistic  and  performance  analysis  of  MSDD. 

•  Performance  comparison  of  OC  and  MSDD. 
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In  Chapter  3,  We  present  a  new  method  to  derive  the  closed-form  expressions  for  the 
exact  BEP  for  OC.  We  first  derive  the  BEP  conditioned  on  the  fading  of  the  interference, 
then  derive  the  unconditional  BEP.  The  expressions  are  for  multiple  interferers  and  multiple 
reception  branches,  with  the  number  of  interferers  less  than  that  of  reception  branches.  The 
modulation  is  BPSK  and  the  channels  are  independent  Rayleigh  fading  channels. 

In  Chapter  4,  We  derive  simpler  asymptotic  expressions  for  BEP  of  OC  for  M-PSK 
modulation  and  one  interferer. 

In  Chapter  5,  we  develop  a  detector  exploiting  MSDD  technique.  The  channel  gain  of 
the  desired  signal  are  assumed  unknown.  M-DPSK  modulation  is  employed.  The  decision 
statistic  for  the  detector  is  derived  based  on  the  principle  of  maximum  likelihood  sequence 
detection  (MLSD).  The  performance  of  the  detector  is  demonstrated  by  simulation  results 
and  analysis  results. 

In  Chapter  6,  another  kind  of  MSDD  is  presented  to  suppress  cochannel  interference 
in  communication  systems.  The  channel  gain  of  the  desired  signal  and  the  channel  phase 
of  the  interference  are  assume  unknown,  but  the  channel  amplitude  of  interference  is  as¬ 
sumed  known  at  the  receiver.  The  interference  signal  is  assumed  to  have  the  same  M-DPSK 
modulation  as  the  desired  signal.  A  maximum  likelihood  sequence  detector  is  developed 
for  detecting  both  the  desired  signal  and  the  interference  signal.  It  is  a  kind  of  multiuser 
detector  employing  MSDD. 

Summary  and  future  work  are  presented  in  Chapter  7. 
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Chapter  3 


Performance  Analysis  for  OC  with 
BPSK  Modulation  and  Multiple 
Interferers 

3.1  Introduction 

As  mentioned  in  Chapter  2,  for  wireless  communication  systems  with  reception  diversity, 
optimum  combining  (OC)  is  a  well-known  approach  to  combat  fading  and  suppress  cochannel 
interference.  It  combines  the  output  of  the  reception  branches  in  an  optimum  way  and 
achieves  the  maximum  output  signal- to-interference-plus-noise  ratio  (SINR). 

Performance  analysis  of  OC  has  been  an  active  research  area.  A  lot  of  effort  has  been 
focused  on  obtaining  closed-form  expressions  for  BEP.  For  the  case  of  only  one  interferer 
existing  in  the  system,  some  work  has  been  done  and  useful  expressions  have  been  derived.  [8] 
and  [14]  gave  out  expressions  for  BEP  when  the  interference  is  not  fading  and  the  modulation 
is  BPSK.  The  non-closed-form  expression  for  fading  interference  was  presented  in  [15],  and 
closed-from  expression  for  BEP  was  obtained  in  [16].  A  non-closed-from  expression  of  BEP 
for  M-PSK  was  shown  in  [17]. 

For  the  case  where  the  number  of  interferers  is  equal  to  or  greater  than  the  number  of 
reception  branches  and  thermal  noise  is  negligible,  BEP  could  be  expressed  in  closed  form 
[18]. 

When  there  are  more  than  one  interferers  and  the  noise  is  not  negligible,  no  closed-from 
expression  for  the  exact  BEP  is  available  until  now,  but  the  performance  has  been  studied 
extensively  through  simulation  [8],  upper  bounds[19],  approximate  expression  [20]  and  non 
closed-from  expressions  (  [21],  [22]).  The  related  work  about  OC  was  summarized  in  [17]. 

The  conventional  way  of  deriving  the  expression  for  BEP  is  to  first  derive  the  probability 
density  function  of  the  SINR  conditioned  on  the  fading  of  the  interference.  The  conditional 
BEP  is  expressed  as  a  function  of  SINR.  The  unconditional  BEP  (which  is  what  we  want) 
is  obtained  by  averaging  the  conditional  BEP  first  over  SINR,  and  then  over  the  fading  of 
the  interference.  Since  it  involves  two  steps  of  averaging,  a  closed-form  expression  for  BEP 
is  difficult  to  obtain. 

In  this  chapter,  we  present  a  new  method  to  derive  BEP.  We  first  derive  the  BEP  con- 


14 


ditioned  on  the  fading  of  the  interference  directly,  then  derive  the  unconditional  BEP  by 
averaging  the  conditional  BEP  over  the  fading  of  the  interference.  This  approach  involves 
only  one  averaging.  Though  complicated,  this  method  provides  closed-form  expression  at 
the  end.  The  basic  configuration  for  our  derivation  is  BPSK  modulation,  Rayleigh  fading 
channel.  We  assume  the  number  of  interferers  is  less  than  the  number  of  reception  branches. 
We  should  point  out  that  the  method  developed  here  could  be  applied  to  the  case  where  the 
number  of  interferers  is  equal  to  or  greater  than  the  number  of  reception  branches. 


3.2  System  Model 

The  system  model  used  in  this  chapter  is  the  same  as  that  mentioned  in  Section  2.2,  which 
is  represented  by  the  following  expression: 

N, 

rk,t  -  JPsCisk  +  ^2  VP'd^k  +  nk,t,  (■  =  1,  •  •  • ,  L  (3-1) 

»=i 

where  the  parameters  are  defined  in  Section  2.2.  We  assume  the  number  of  independent  re¬ 
ception  branches  L  is  greater  than  the  number  of  interferers  AT/.  sk  is  the  desired  transmitted 
BPSK  symbol.  The  interference  signal  shk  is  assumed  to  be  Gaussian  distributed,  with  zero 
mean  and  unit  variance. 

The  channel  gains  ct  of  the  desired  signal  and  of  the  interference  are  assumed  to  be 
independently  and  identically  distributed  (i.i.d.),  zero-mean,  circularly  symmetric,  complex 
Gaussian  random  variables  (Rayleigh  fading),  with  variance  1/2  per  dimension. 

As  mentioned  in  Section  2.2,  the  signal  model  in  vector  notation  is 

N, 

rk  =  y/KcSk  +  v/f^X>,fc  +  n*.  (3.2) 

1=1 

3.3  Derivation  of  Conditional  BEP 

As  shown  in  Fig.  2.4,  for  OC  detector,  the  received  signal  rk  are  weighted  and  combined  to 
get  an  output  signal.  The  weighted  vector  that  yields  the  largest  SINR  is  [8] 

w  =  R^c.  (3-3) 

The  output  of  the  combiner  is  wHrk.  For  BPSK  modulation,  the  decision  rule  of  the  detector 
is:  if  Re(wffrfc)  ^  0,  the  decision  is  made  that  1  is  transmitted;  otherwise  the  decision  is 
made  that  —1  is  transmitted. 

From  now  on,  we  assume  the  transmitted  symbol  is  sk  =  1,  then  from  (3.2),  the  received 
signal  is 

_  Nl 

Tk  =  s/PsC  +\/PlJ2  CiSi’k  +  (3‘4) 

t=l 
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Define  random  variable 


D  -  2 Re  (w^r*)  =  wHrk  +  (wHrk)* .  (3.5) 

Note  that  D  is  random  due  to  that  channels  c  and  c /,  interference  Si<k  and  noise  n*  are 
random.  According  to  the  decision  rule,  when  D  <  0,  the  decision  is  made  that  —1  is 
transmitted  and  an  error  occurs.  Therefore  the  BEP  is  Pe  =  Pr(D  <  0).  Our  purpose  is  to 
expression  Pe  as  a  closed-from  function  of  the  SNR,  SIR,  number  of  reception  branches  L 
and  number  of  interferers  Nj. 

It  is  very  difficult  to  derive  the  expression  for  Pe  directly.  We  derive  it  in  two  steps.  In 
this  section,  we  treat  the  fading  of  the  interference  c /  as  constant  (which  is  actually  random). 
More  specifically,  in  this  section,  c /  =  [ci,  C2,  •  •  •  ,  c^]  is  considered  as  non-random,  rk  (and 
consequently  D )  is  random  only  due  to  the  randomness  of  c,  si)fc  and  nfc.  Based  on  this 
assumption,  we  derive  the  conditional  BEP  P  (e|c/)  =  Pr(£)  <  0|c/).  In  the  next  section, 
we  will  consider  the  random  effect  of  C/  and  derive  the  unconditional  BEP  Pe  by  averaging 
the  conditional  BEP  P  (e|cj)  over  c /. 

When  the  fading  of  the  interference  C/  is  treated  as  constant,  the  covariance  matrix  of 
the  interference  plus  noise  zk  is 


Ni 

Hz  =  Pi  ^  Cjcf  +  o2Il  ,  (3.6) 

i=l 

where  the  superscript  H  denotes  the  Hermitian  transposition,  a 2  is  the  power  of  the  noise 
and  1 1  is  an  identity  matrix  of  rank  L. 

Let  *&D(ju)  be  the  characteristic  function  of  D  conditioned  on  c j.  Using  similar  procedure 
as  in  [7,  Appendix  B],  it  can  be  shown  that  the  conditional  BEP  is 

P(e|c/)  =  Pr(£)  <  0|cj) 

=  __L  r*“  Wi, 

2 *3  J-oo+je  W 

—  —  ^2  R-es 

lm(wm)>0 

where  e  is  a  small  positive  number;  Res  ^ denotes  the  residue  of  function 
at  pole  um.  The  summation  is  taken  over  the  poles  in  the  upper  half  of  the  complex  plane. 

The  task  at  hand  is  to  evaluate  the  characteristic  function  $d(jw).  To  that  end,  we  try 
to  express  D  in  (3.5)  as  a  quadratic  form  of  complex-valued  Gaussian  random  variables  and 
apply  the  results  in  [7,  Appendix  B]. 

Substitute  (3.3)  into  (3.5), 

D  =  c^r*  +  (c*R ?rky  .  (3.8) 

Diagonalize  Rz  as  R*  =  UzAzUf,  where  Az  =  diag  (At,--  -  ,AL),  Aj,---,Al  are  the 
eigenvalues  of  Rz,  and  Uz  is  the  unitary  matrix  whose  columns  are  the  eigenvectors  of  Rz. 
According  to  [22],  we  can  assume  that  the  L  eigenvalues  of  Rz  satisfy:  Ai  ^  A2  ^  ^ 


$d(. ju) 


UJ 


(3.7) 
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^Ni  v2  and  Am  —  o 2  for  m  =  Nj  +  1,  A/j  +  2,  •  •  •  ,  L.  Define  the  vector  of  the  non-zero 
eigenvalues  as  A  =  (Ai,  •  •  •  ,  Ajy>) . 

From  Rz  =  UZAZU^  we  have 


r:1  =  U-A^U; 


-IttH 


Substitute  the  above  expression  into  (3.8), 


D  = 


c"UA:'V"rk  +  (c"U,A;1Ufrt)‘ 

(U?c)"a;>  (Ufr.)  +  [(Ufc)"  A-  (Ufr*)]'. 


(3.9) 


(3.10) 


Define  the  whitened  interference-plus-noise  vector  x  =  [xi,X2,-“  ,xL)T  and  modified 
channel  vector  g  =  [gi,g2,  •  •  •  ,9l}T  as 


x  = 


g  = 


Note  that  since  U2  is  unitary,  the  vectors  g  and  c  have  the  same  distribution. 
Substitute  (3.11)  and  (3.12)  in  (3.10), 


D  =  g^x+Cg^x)’ 

L 

=  }  y  \n  {9mXm  +  9mxm)  • 

m= 1 


(3-11) 

(3.12) 


(3.13) 


In  (3.13),  D  is  a  quadratic  form  of  complex-valued  random  vectors  x  and  g.  Applying  the 
results  in  [7,  Appendix  B]  to  (3.13)  (see  Appendix  A  in  this  report  for  details),  we  obtain 
the  characteristic  function  of  D  as 


2  L 

$dUu)  =  n  i 

m=l 


i _ 


(3.14) 


where 


a?  +  \/ft  +  i 

m=  1, 

-”,Ni 

l^m  —  ^ 

/_/»  +  i 

^m-Nj  y  \n-Nj  Am_ 

Ni 

m  =  Ni  + 

!»•••  ,2A/> 

(3.15) 

/J00 

m  — 

2Ni  +  1,  •  •  • 

,2  NT  +  (L- 

-AT,) 

k  H  01 

m  = 

:  2Ni  +  (L  — 

Ni)  + 

,21 

"ft 

o 

o 

ll> 

VE, 

a 2 

(3.16) 

•ft 

o 

ll> 

VPs 

<r2 

(3.17) 
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Notice  that  with  regard  to  the  fading  of  the  interference  c j,  noo  and  /i0i  are  constants.  They 
only  change  with  Ps  and  a2.  Define  the  signal  to  noise  ratio  (SNR)  7  as 


Ps 

7=^- 

(3.18) 

Then  p0o  and  poi  in  (3.16)  and  (3.17)  could  be  written  as 

Poo  -  ^5- (7+VV  +  7) 

(3.19) 

fo,  -  ^(7  V¥T^). 

(3.20) 

Using  (3.15)  and  (3.14),  we  have 


•Mj'w) 


uj 


1 

UJ 


’2  Nt  . 

rr _ I _ 

i1  - 


(1  -  jtioou)L  Nl  (1  -  jpoiw) 


L—Nj  • 


(3.21) 


Since  nm  <  0  (m  =  Nj  +  1,  Nr  +  2,  •  •  •  ,  2iVj)  and  /x0i  <  0,  for  $D(joj)/uj  in  (3.21),  the 
poles  -j/nm  (m  =  Nj  +  1,  N/  +  2,  •  •  •  ,  2 TV/)  and  pole  -j/nm  are  in  the  upper  half  of  the 
complex  plane.  Therefore  (3.7)  becomes 


2  N{ 


P  (e|c/)  =  —  X  Res 

*71=^/  +  ! 


^z?(jw)  1 

j  J 

—  Res 

$dC?w)  1 

;  3 

UJ  Hm . 

u  //01 

(3.22) 


Substituting  (3.21)  into  (3.22),  and  carrying  out  the  calculation  of  the  residues  (Appendix 
B),  we  get  the  conditional  BEP  as 


P(e  |cf) 

2Ni 

=  -  E 

m=Nj+ 1 

1 


(  m— 1 


u2L 

Fm 


- 


n 


1 


,n=l  V™/  \n=m+l/Xn 

1 


(POO  -  Mm)1”  ^  (Pol  -  Hm)L  Nl\ 

_ ! _ 1 _ y' 

(mooMoi)L~W/  (L  —  Ni  -  1)!  '  £ 


L-N/-1 


(L  —  Ni  +  £  -  1)! 

P 


2ATj 

+E 


(A^n/^Ol) 


L-Nj-i 


(Pn  -  M0l)L 

(PooPoi)l_JV/+< 


(POO  —  P01) 


(3.23) 


P  (e|c/)  is  a  function  of  fim  (m  =  1, 2,  •  •  •  ,  2Nj).  Since  pm’s  (m  =  1, 2,  ■  •  •  ,  2JV/)  are  functions 
of  eigenvalues  Am’s  (m  =  1, 2,  •  •  •  ,  A/),  P  (e|cj)  is  a  function  of  A  =  (Ai,  •  •  •  ,  A Nj) ,  which  is 
related  to  the  fading  of  interference  Cj. 
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Special  case:  no  interference 

We  can  use  the  method  detailed  above  to  derive  the  BEP  for  the  case  where  there  is  no 
interference.  In  this  case,  OC  becomes  maximum  ratio  combining.  Since  the  power  of 
interference  Pj  =  0,  from  (3.6)  we  have  the  covariance  matrix  Rz  =  cr2I l-  The  eigenvalues 
of  Rz  would  be  Am  =  a2  for  m  =  1, 2,  •  •  •  ,  L.  Referring  to  (3.15),  (3.16)  and  (3.17),  we  have 
/zm  =  /zoo  for  m  =  1, 2,  •  •  •  ,  Nj  ;  Hm  =  /toi  for  m  =  AT/  +  1,  Nj  +  2,  •  •  •  ,  2A/j.  (3.21)  becomes 


$z>(jw)  _  1 


w  to  (1  -  jfiooto)1  (1  -  jno\u)L 
Substitute  the  above  expression  in  (3.7),  we  get  the  BEP  as 


Pe  =  —Res 


.  1 

to  ’  J  Hoi. 


(3.24) 


(3.25) 


Since  there  is  no  interference,  the  BEP  is  unconditional  BEP.  After  carrying  out  the  calcu¬ 
lation  of  the  residues,  we  get 

VL-l 


P,= 


(—Voi) 


(Hoo  ~  Hoi )L{L  ~  1)! 


E 

,z=o 


(L-l+^)! 


Moo 


&  (Moo  —  Moi)* 


(3.26) 


(3.26)  is  the  expression  of  BEP  for  without  interference.  We  can  express  it  in  a  more 
convenient  form. 

Substitute  (3.19)  and  (3.20)  in  (3.26),  we  have 


1  1 
(L  —  1)! 


(L-l+^)! 


(3.27) 


Define 


Then  (3.27)  becomes 
P, 


ii 


(1  -  h) 


gmiiM')' 


(3.28) 


(3.29) 


which  is  the  same  as  Eq.  (14-4-15)  in  [7].  We  derive  (3.29)  without  the  need  of  integration. 
In  [7],  (3.29)  was  obtained  by  integration. 

When  SNR  7  >  1,  /z  =  (l  +  2«1-^.  (3.29)  could  be  approximated  as 


3.4  Derivation  of  Unconditional  BEP 


The  unconditional  BEP  Pe  is  obtained  by  averaging  the  conditional  BEP  P  (e|cj)  over  the 
fading  of  the  interference  c j.  Since  P  (e|cj)  is  a  function  of  the  vector  of  eigenvalues  A,  and 
A  depends  on  c /,  averaging  P(e|c/)  over  cj  is  equivalent  to  averaging  it  over  A.  Hence 


pe  =  J  P(e\ci)px(X)dX 

where  px(y>  is  the  joint  probability  density  function  of  A  given  out  in  [22]  as 


PxW  =  Kojp 


« (^r 


X 


(3.31) 


(3.32) 


where 


K0 


1 

n, "',(£-*)!  nsw-i)!' 


(3.33) 


In  (3.23),  the  conditional  BEP  P(e|c/)  is  a  function  of  /im’s.  As  shown  in  (3.15),  pm' s 
are  irrational  functions  of  eigenvalues  Am’s.  Since  integration  of  irrational  function  is  very 
difficult  to  carry  out,  we  try  to  convert  the  conditional  BEP  P  (e|cj)  into  a  rational  function 
of  some  other  variables  and  then  averaging  it  over  those  variables.  For  this  purpose,  we 
define  constant  (with  respect  to  c /)  rj  as 


T]  = 


(3.34) 


and  variables  ym  as 

ym  =  J^+l,  m  =  l, 2,  —  ,N,.  (3.35) 

Also  define  vector  y  =  [j/i,  j/2>  •  •  •  ,  J/n>]  •  Since  A  =  (Ai,  •  •  •  ,  A  at,)  is  a  random  vector,  so  is  y. 

From  (3.34)  and  (3.35)  we  have 

a2  =  Ps  (rf  —  l)  (3.36) 

Am  =  P«  (l/m  -  l)  m  =  l,2,“-  ,Nj.  (3.37) 

Using  (3.34),  (3.16)  and  (3.17),  we  get 

Poo  =  — y  +  W 

a* 

VPs  n  v 
Poi  =  (!“»?)■ 
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(3.38) 

(3.39) 


From  (3.35)  and  (3.15),  we  have 


Mm 

Mm+Nj 


y/Ps 


(l  ?/m) 
(1  —  Vm) 


(3.40) 

(3.41) 


for  ra  =  1, 2,  •  •  •  ,  Nj. 

Substituting  (3.38),  (3.39),  (3.40)  and  (3.41)  into  (3.23),  after  some  straightforward  ma¬ 
nipulations,  we  get  the  conditional  BEP  as  a  function  of  the  variables  ym’s, 


^(e|y) 


Ni 

-  53  /m  (y) 

m—  1 
L-AT/-1 

E 

e=o 


{L  —  Nj  —  1)! 


(-1) 


L-AT/-1 


(L  —  Ni  +  £—  1)! 


AT/ 

i+E^ 

m= 1 


L-N!-e 


(2*7) 


L-Ni+e 


H)‘ 


(3.42) 


where 


/m(y) 


i-ym  (i~t) 

2fc  te-v2)1""' 


and 


5m,/  (y)  =  (-1) 

Fl-nt 


L—Ni—t  (1  +  r?)L 


”  2 Vm  (yg,  - 

'“{BrSMAr*} 


with 


FL-Nl-t  (ym)  =  -  (1  +  2/m)  (f?  -  Vm)1  Nl  1  +  (1  -  5m)  (*?  +  2/m)L  ^ 


(3.44) 


(3.45) 


Obviously,  conditional  BEP  P  (e|y)  is  a  rational  function  of  the  random  variables  ym's. 
The  joint  probability  density  function  of  y  is 


JVty)  =  „»( \)p--pL 

VXK  I  dyx  dyN[ 


A=y 


X 


,2\2 


n  w-p2) 

Ll<i<7<W/-l 


y\V2  ■  •  •  yjv/ 


(3.46) 
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for  yi  ^  y2  ^  ^  Vn,  >  y,  where 


0  =  Ps/Pi 


is  the  signal  to  interference  ratio  (SIR)  and 


K  = _ f_ _ pur, 

[n£,  (i  -  *•)!]  [n£,  w  -  ••)!] 


(3.47) 


(3.48) 


The  unconditional  BEP  Pe  is  obtained  by  averaging  the  conditional  P  (e|y)  over  the 
random  vector  y, 


p,  =  J  P(e|y)py(y)<iy 

=  -  E  //» (y)jv(y)*r  -  (L-^,_i)i  <_1> 

m=l 


L-Nj-l 


L—Nj—l 


h  v  v  V 


i  +  9m^(y)py{y)dy 


L—Nj—t 


-fy+t  * 


(3.49) 


In  order  to  carry  out  the  averaging  of  gm<t  (y)  >  we  need  to  express  it  in  a  more  convenient 
way.  Specifically,  for  easy  to  integrate,  we  express  the  function  (t/m)  in  9m, e  (y) 

(shown  in  (3.44)  and  (3.45))  as  a  sum  of  (y^  -  y2)  to  integer  power.  From  Appendix  C,  we 
have 


[(£—*>—/)/ 2] 


Fl-Ni-1  (Vm)  —  2ym  ^  aL-N[-t,t  (ym  y  )  j 


(3.50) 


where  [(L  —  iV/  —  ^)  /2]  denotes  the  largest  integer  that  is  equal  to  or  less  than  (L  —  Ni  —  i)  /2, 
and  is  evaluated  as: 


a L-Ni-t 


\(L  —  Nj  —  i  —  \  — 1\  .  (L-Ni-£-l-t\] 

-  [{  t  )(1-’,)-2’,(  t-i  )\ 


\L-Ni-l~l-2t 


(3.51) 


When  calculating  aL-Nl-tyU  we  assume  (™)  =  0  for  m  <  n  or  n  <  0. 

The  first  sum  of  multiple-fold  integrals  Y,m= i  /  /"*  MPyM^y  in  (3.49)  is  evaluated  in 
Appendix  D.  The  final  results  is: 


fV/  r  ,  Ni-lNi-l 

E  /  /-MMyMy  =  (--)  E  E  (-D*  »(!>.«) 

m-l  ^  '  ''  p=0  9=0 


(3.52) 
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where  Bq  (5  =  0, 1, 2,  •  •  • )  is  a  series  that  could  be  evaluated  as  follows: 


^exp  (£r?2)  Q  (v/2/ty) 

(3.53) 

i-1)2  exp  (jhf)  Q(x/2f)ri'j 

(3.54) 

[(29-1) 

w 

-if  B,-i  + 

(3.55) 

and  when  Nr  =  1,  H  (p,q)  =  (X^ryyr;  when  Nj  >  1, 

H(p,q ) 

_ 1 _ 

”  [n£.  (i  -  «•)!]  [n£,  (N<  -  os] 

e  e  detw-  <3-56) 

mi+m2+-4mL>i=L-l-p  ni+n2H - Vni-i=L-\-q 

mi€{0,l}  ni€{0,l} 


where  det  W  is  the  determinant  of  a  (Nj  —  1)  x  (jVj  —  1)  matrix  whose  element  on  the  «-th 
row,  j-th  column  is 


Wij  =  ( rrij  +  nj  +  L-NI  +  i+j-  2)!. 


(3.57) 


The  summations  in  (3.56)  are  taken  over  all  sets  of  indices  meeting  the  stated  conditions. 
H  (p,  q)  depends  on  L,  Ni,  p  and  q,  and  is  independent  of  SNR  7  and  SIR  /?.  Therefore  we 
can  calculate  it  for  once  and  then  use  it  for  all  value  of  SNR  and  SIR. 

Similarly  as  in  Appendix  D,  we  can  get  the  second  sum  of  multiple-fold  integrals  in  (3.49) 
as: 


pL-Nj+l 

[(L-Nr-t)/2] 

aL~Ni~t’t  pe+t+i 

Nj-l  Nf-l  /  ,o\  p 

EE  (-iyH{p,q)(t  +  t  +  q)\l^)  .  (3.58) 

p= 0  q- 0  \1/ 


(3.49)  is  the  expression  for  the  exact  BEP  of  OC  for  BPSK  over  Rayleigh  fading  channels. 
Prom  it  and  other  related  expressions  we  can  calculate  BEP  for  any  given  number  of  branches 
L  ,  number  of  interferers  Nj  (with  Nj  <  L),  SNR  7  =  Ps/cr2  and  SIR  /?  =  Pa/P/. 

Based  on  (3.49),  we  can  derive  simpler  asymptotic  expression  for  the  special  case  SNR 
7  »  1  and  SIR  /5<1. 
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Special  case:  SIR  (3  <c  1  and  SNR  7  1 

Notice  that  (3.52)  and  (3.58)  all  contain  the  term  pL~Ni+l.  Since  0  <  1  and  L  >  NJt  f}L~N'+l 
0.  We  conclude  that 


Y  [  fm(y)py(y)dy  ~  0 

m=l  “ 

N,  r 

Y  /  9m,i(y)py{y)dy  »  0. 

m=l  " 


(3.59) 

(3.60) 


Prom  (3.49)  we  have 
1 


e~  (L  —  Ni  —  1)!  j  t? 

Since  7  >  l,from  (3.34), 


L-N,-  1 


(L  —  TV/  +  ^  —  1) 


!H)‘ 


L-Ni-t 


(2*7) 


L-AT/+< 


-.  (3.61) 


/I  1 

r)=\l-  +  1~1+  2^ 


Substitute  (3.62)  in  (3.61), 


(L-TV/-1)! 


(-1) 


L-Nt-I 


L-Ni-1 

E 

/=o 


(L-iVj +  /-!)! 


!B)‘ 

L-AT/-1 

E 


(-*) 


L-N!-l 


,  L-Nt+l 


(2) 


(L  —  Ni  + 1  —  1)! 


-  ^l-NijL-Ni  (L  —  Ni  —  1)!  ^ 

1  /2  (L  —  TV/)  —  1\ 

-  ^L—NjyL—Nj  ^  (L- NT)  )’ 


(3.62) 


(3.63) 


Comparing  (3.63)  with  (3.30),  we  can  see  that  the  performance  of  a  system  with  L 
branches  and  JV>  large  interferers  is  equivalent  to  that  of  a  system  with  (L- Nr)  branches 
but  without  interference. 


3.5  Numerical  Results 

As  mentioned  in  the  introduction  section,  the  performance  of  OC  has  been  evaluated  in 
several  papers  (without  using  closed-form  expressions).  We  don’t  want  to  repeat  that.  We 
only  use  numerical  results  to  show  that  the  analysis  results  from  the  closed-form  expressions 
match  the  simulation  results.. 

Fig.  3.1  to  Fig.  3.4  show  BEP  versus  SNR  for  different  SIR.  Fig.  3.1,  3.2  and  3.3  are  for 
L  =  4  branches,  and  TV/  =  1, 2, 3  interferers,  respectively.  Fig.  3.4  is  for  L  =  8  and  Nr  =  5. 

In  all  the  figures,  the  analysis  results  are  very  close  to  the  simulation  results.  That  proves 
the  validity  of  the  analytical  expression  for  BEP. 
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Bit  Error  Probability 


V  SIR=0.02  (Simulation) 

-  SIR=0.02  (Analysis) 

A  $IR=1  (Simulation) 

-  -  SIR=1  (Analysis) 

$  StR=10  (Simulation) 
••••  SIR=10  (Analysis) 

*  SIR=50  (Simulation) 

•  -  •  SIR=50  (Analysis) 


Figure  3.1:  BEP  versus  SNR  for  L  =  4  branches,  iVj  =  1  interferer. 


Figure  3.2:  BEP  versus  SNR  for  L  =  4  branches,  Nr  =  2  interferers 


Bit  Error  Probability 


SNR  (dB) 

Figure  3.3:  BEP  versus  SNR  for  L  =  4  branches,  JV/  =  3  interferers. 


SNR  (dB) 

Figure  3.4:  BEP  versus  SNR  for  L  =  8  branches,  Ni  =  5  interferers. 
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Chapter  4 


Asymptotic  BEP  for  OC  with  M-PSK 
Modulation  and  One  Interferer 


4.1  Introduction 

In  Chapter  3,  we  derived  closed-form  expressions  of  BEP  for  OC.  The  expressions  were 
for  BPSK  modulation,  Rayleigh  fading  channels,  multiple  interferers  with  the  number  of 
interferers  less  than  the  number  of  reception  branches. 

An  expression  was  given  out  for  the  BEP  with  M-PSK  modulation  and  one  interferer  in 
[17].  It  involved  integration.  In  this  chapter,  we  will  derive  closed-form  expressions  of  BEP 
for  asymptotically  high  SNR.  The  asymptotic  expressions  give  intuitive  inside  of  OC  and  are 
easy  to  calculate.  We  need  these  expressions  to  compare  the  performance  of  OC  with  that 
of  MSDD  later. 

Expressions  were  given  out  for  the  BEP  of  OC  with  M-PSK  and  multiple  interferers  in 
[22].  Future  work  needs  to  be  done  to  get  asymptotic  expressions  from  that. 


4.2  Existing  Expressions 


For  M-PSK  signals,  the  conditional  symbol  error  probability  (SEP)  for  OC  is  [17,  Eq.  10.29] 
(some  notations  have  been  changed  to  agree  with  this  report) 


s, M-PSK 


W.)  =  \  [ 

ft  J  0 


(M— 1  )ir/M  r 


^7f|Ai  (5)|5_  sin2(jr/Af) 


sin2  6 


dd, 


(4.1) 


where  function  M7<  (*,(«)  is  the  moment  generating  function  of  SINR  jt  conditioned  on  Ai, 
the  largest  eigenvalue  of  the  interference  plus  noise  covariance  matrix  Rz,  and  integer  M  is 
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the  number  of  symbols  in  M-PSK  modulation.  Prom  [17,  Eq.  10.11],  we  have 

(s)|s=_«in£MMl 

sin2  $ 

sin2  (7 r/M)^ 
sin2 

sin2  # 


=  M. 


(sii 

sin2  0  ) 


f  ( 


sin2  0 


sin2  9  +  sin2  (n/M)  7  /  \^sin2  9  +  sin2  (7 r/M)  &  J 


(4.2) 


The  unconditional  SEP  Ps, m-psk  is  obtained  by  averaging  the  conditional  SEP  over  the 
random  eigenvalue  Ax, 

P«.m- psk  =  /  Ps.m-psk(P|Ai)??Ai  (Ai)dAi,  (4.3) 


poo 

Ps, M-PSK  =  /  Ps,M-PSK(P|Ai)pAl(Ai)dAi, 

J  <r2 


where  pAl(Ax)  is  the  probability  density  function  of  Ai.  Prom  (3.32)  in  Chapter  3,  for  where 
there  is  only  one  interferer  (i.e.,  AT/  =  1),  we  have 


/x  x  11  (  Ax  -a2\  ( Ai  -<t2n\L  1 

PAl(Al)_  (L-l)!P/eXP  V  Pi  A  Pi  ) 


(4.4) 


For  Gray  coding  and  large  SNR,  the  relation  between  SEP  Ps, m-psk  and  BEP  Pb> m-psk  is 

1  „ 


Pb, M-PSK 

Prom  (4.5),  (4.3)  and  (4.1),  the  BEP  is 


log 2M 


-s, M-PSK- 


(4.5) 


Pb, M-PSK  « 


<4-6) 

(4.6)  is  the  expression  for  the  general  cases.  The  approximation  is  due  to  the  transfor¬ 
mation  from  SEP  to  BEP. 

Next  we  derive  the  asymptotic  BEP  of  M-PSK  for  high  SNR  7  1,  for  the  cases  of  no 

interference  and  SIR  -C  1. 


4.3  Asymptotic  Expression  for  No  Interference 

When  there  is  no  interference,  P/  =  0.  The  OC  is  equivalent  to  maximum  ratio  combining 
(MRC).  In  this  case,  Ai  is  a  constant,  Ax  =  <r2.  Substitute  it  into  (4.2), 

(4.7) 


/  sin2  _  ( _ sin2# _ \ 

7^Al  V  sin2#  J  \sin2 #  +  sin2 (tt/M) 7 / 


(4.6)  becomes 


(M— 


M. 


7t|Al 


(- 


sin2  (tt/M)\ 
sin2  #  / 


dO. 


(4.8) 
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sin2  (7 r/M)\ 


For  high  SNR  7  1, 

^  sin2  6  )  7Lsin2i  (n/M) 

Substitute  (4.9)  into  (4.8), 


sin2i0. 


1  1  1  f 

g 2M  7r  7l  sin2t  (ir/M)  Jo 


sin2L  Odd. 


b,M-PSK  ~  \0g2M  7r-y^sin2t  (ir/M) 

The  above  integration  is  treated  separately  for  M  =  2  and  M  ^  4. 

4.3.1  BPSK,  M  =  2 

Using  Eq.  (2.513)  in  [23],  we  have 


if 


*/2  1  /2L 

sin2L  Odd  = 


22L+1 


(?) 


Since 


O- 


(2L)l 
LILl 
2L  (2L  —  1)! 
“  L  (L  -  1  )\L\ 

-  ■(?:.')• 


ir™2Lm=H2L-i)- 


Substituting  (4.13)  into  (4.10),  we  get  the  BEP  for  BPSK  as 

■Pb.BPSK  ~  2 1 


,/2L-l\  1 

'  \  L  -  1  )  4l7l  ’ 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


(4.13) 


(4.14) 


4.3.2  M-PSK,  AO  4 

When  M  =  4  (QPSK),  define  a  function  p  (L)  such  that 

'  dn“ft»  =  p(L)JE(“_11). 

The  relation  between  p  (L)  and  L  is  shown  in  Fig.  4.1  and  the  following  table. 


(4.15) 


L 

2 

3 

4 

5 

6 

10 

20 

80 

P(L) 

1.92 

1.97 

1.99 

1.99 

2.00 

2.00 

2.00 

2.00 
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Substitute  (4.20)  into  (4.10),  we  get  the  BEP  for  M  ^  4  as 


Pb,M-PSK  ~ 


log2M  jL  sin2L  (jr/M)  22L  \  L  —  1  / 
2  /2 L  -  1\  1 

log2M  \  L  —  1  /  4i71,  sin2L  (ir/M) 


(4.21) 


In  the  numerical  section,  it  is  proved  that  for  high  SNR  7  1,  the  approximate  BEP 

expressions  in  (4.14)  and  (4.21)  are  very  close  to  the  BEP  yielded  by  (4.8). 


4.4  Asymptotic  Expressions  for  SIR  <C  1 


In  [17,  Chapter  10],  it  is  shown  that  the  unconditional  SEP  in  (4.6)  can  be  very  closely 
approximated  as 


■Ps.M-PSK  ~  -Ps,M-PSk(-E'|Ai)|Ai=ji 


1 

7T  Jo 


(  sin2(7r/M)\ 

\  sin2  9  ) 


d9, 


(4.22) 


where  Ai  is  the  expectation  value  of  Ai,Ai  =  LP[  4-  a2.  Substitute  Ai  into  M~,t |Ai  (*)  (from 
(4-2)), 

/  sin2(7r/M)\  _  / _ sin2 9 _ \L  1  / _ s*1*2# _ \ 

7t|Xl  V  sin2  9  )  ~~  \sin2  0  +  sin2  (i r/M)  7  J  ^  sin2  0  +  sin2  (7 r/M)  LPf^  J 

(4.23) 


For  SIR  <  1,  Ps  <C  Pi,  hence  Ps  <  LPi  +  a2.  Applying  this  to  (4.23),  we  obtain 


M. 


7t|*l 


(- 


sin 


(* m\ 

sin2  9  J 


7L-1sin2^i  ^  (tt/M) 


sin 


PL~1'>9. 


(4.24) 


The  difference  between  the  MGF  M7t|Al(-)  for  the  case  of  SIR  <C  1  in  (4.24)  and  the 
MGF  without  interference  in  (4.9)  is  the  loss  of  a  diversity  degree  of  freedom  in  (4.24).  We 
conclude  that  the  BEP  for  this  case  is  obtained  by  replacing  L  in  (4.14)  and  (4.21)  with 
L  —  1.  Therefore  for  optimum  combining  in  the  presence  of  an  interference  source,  when 
SNR  »  1  and  SIR  C  1, 


A.bpsk  ~ 


/2  (L  —  1)  —  l\  1 

\  (L  —  1)  —  1  /  4t_17L_1 


for  BPSK 


(4.25) 


and 


2  /2  (L-  1)  -  1\ _ 1 _ 

log2M  \  (L  —  1)  —  1  /  4x,-17£,_1  sin2^-1^  (n/M) 


for  M-PSK. 


(4.26) 
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Figure  4.2:  Comparison  of  asymptotic  results  and  exact  results,  no  interference,  L  =  4. 

4.5  Numerical  Results 

In  this  section  we  provide  some  numerical  results  to  show  how  close  the  asymptotic  expres¬ 
sions  are  to  the  non-asymptotic  expressions  (which  we  call  ’exact’  expressions). 

Since  the  BEP  for  QPSK  is  very  close  to  that  of  BPSK,  we  only  show  BEP  for  BPSK 
and  8-PSK.  All  the  results  are  for  L  —  4  reception  branches. 

Fig.  (4.2)  show  that  when  there  is  no  interference,  the  asymptotic  results  yielded  by 
(4.14)  (for  BPSK)  and  (4.21)  (for  8-PSK)  are  very  close  to  the  exact  results  yielded  by  (4.8) 
for  SNR  >  15  dB. 

Fig.  4.3  and  (4.4)  show  the  results  for  the  case  with  interference.  The  asymptotic  results 
are  yielded  by  (4.25)  (for  BPSK)  and  (4.26)  (for  8-PSK).  The  ’exact’  results  are  yielded  by 
(4.22)  1.  For  SIR  =  0  dB  and  BPSK  in  Fig.  4.3,  the  asymptotic  results  are  not  very  close  to 
the  exact  results  since  SIR  is  not  much  less  than  1.  For  all  other  cases  shown,  the  asymptotic 
results  are  very  good  approximation  to  exact  results  for  SNR  >15  dB. 


1The  difference  between  (4.22)  and  (4.6)  is  so  small  that  it  could  be  neglected.  Hence  in  here  we  call 
results  from  (4.22)  exact  results. 
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Figure  4.3:  Comparison  of  asymptotic  results  and  exact  results,  L  =  4  branches,  SIR 
dB. 
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Figure  4.4:  Comparison  of  asymptotic  results  and  exact  results,  L  =  4  branches,  SIR 


Chapter  5 

Multiple-Symbol  Differential 
Detection  with  Known  Covariance 
Matrix  of  Interference  Plus  Noise 


5.1  Introduction 

In  the  previous  two  chapters,  we  analyzed  the  performance  of  OC.  OC  is  a  coherent  detection 
technique  which  requires  the  channel  phase  of  the  desired  signal.  In  this  and  the  next  chap¬ 
ter,  we  present  a  noncoherent  detection  scheme  which  is  called  multiple  symbol  differential 
detection  (MSDD). 

MSDD  was  first  proposed  for  detecting  M-PSK  signals  transmitted  over  an  additive 
white  Gaussian  noise  (AWGN)  channel  [24].  The  main  advantage  of  MSDD  is  that  it  does 
not  require  a  coherent  phase  reference  at  the  receiver  (it  does  require  however,  the  ability  to 
measure  relative  phase  difference). 

MSDD  performs  maximum  likelihood  detection  of  a  sequence  of  information  symbols 
based  on  a  finite  observation  interval.  The  method  was  presented  as  a  bridge  of  the  gap 
between  the  performance  of  coherent  detection  of  M-PSK  and  conventional  differential  de¬ 
tection  of  M-ary  differential  phase  shift  keying  (M-DPSK).  The  channel  phase  was  assumed 
to  be  constant  over  multiple  symbol  intervals  and  was  unknown  to  the  receiver.  In  [24]  it 
was  shown  that  for  a  long  observation  interval,  the  performance  (in  terms  of  the  required 
SNR  for  a  given  BEP)  of  MSDD  approached  that  of  coherent  detection  (with  differential 
encoding  at  the  transmitter). 

MSDD  was  extended  to  trellis  coded  M-PSK  in  [25].  MSDD  for  the  fading  channel  was 
analyzed  in  [26]  and  for  correlated  fading  in  [27].  MSDD  application  to  multiuser  CDMA 
was  considered  in  [28].  Performance  of  MSDD  with  narrow-band  interference  over  nonfading 
channel  was  discussed  in  [29].  A  system  with  MSDD  and  reception  diversity  was  formulated 
in  [30]  and  [31],  while  [32]  considered  MSDD  with  transmit  diversity. 

In  this  chapter,  we  derive  an  extension  to  MSDD  for  communication  in  the  presence  of  a 
single  interference  source.  The  channel  of  the  desired  signal  is  a  diversity  Rayleigh  channel 
with  multiple  outputs.  The  channel  realizations  at  each  output  are  mutually  independent, 
constant  over  the  observation  interval  and  unknown  to  the  receiver.  The  Gaussian  assump- 
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tion  is  made  with  respect  to  the  aggregate  of  interference  plus  noise.  The  covariance  matrix 
of  the  interference  plus  noise  is  assumed  known.  The  MSDD  decision  statistic  is  derived 
based  on  the  principle  of  maximum  likelihood  sequence  detection  (MLSD).  A  closed-form 
expression  for  the  pairwise  error  probability  (PEP)  is  derived.  A  closed-form  expression  for 
the  BEP  is  intractable,  however  one  is  obtained  for  an  approximation  to  the  union  bound. 
The  approximation  utilizes  only  dominant  terms  in  the  union  bound  and  it  is  shown  to  be 
a  good  approximation  of  the  BEP.  We  show  that  with  an  increasing  number  of  symbols  in 
the  observation  interval,  the  performance  of  MSDD  approaches  that  of  OC  (with  differential 
encoding  at  the  transmitter). 

In  the  course  of  designing  simulations  for  evaluating  MSDD,  we  realized  that  there  was 
no  efficient  MSDD  algorithm  available  for  MSDD  with  diversity.  The  computational  com¬ 
plexity  of  direct  computation  of  the  decision  statistic  grows  exponentially  with  the  number 
of  symbols  in  the  observation  interval.  For  single  channel  MSDD,  an  optimum  algorithm 
was  proposed  in  [33].  Sub-optimal  decision  feedback  algorithms  for  the  single  channel  case 
were  suggested  in  [34],  [35],  and  [36].  In  this  chapter,  we  modify  the  sub-optimal  decision 
feedback  algorithm  in  [36]  for  application  to  MSDD  with  diversity.  The  main  improvement 
over  published  algorithms  is  the  introduction  of  iterations  for  symbol  detection. 


5.2  System  Model 

For  convenience,  we  repeat  the  system  model  which  was  initially  described  in  Section  2.2, 
since  there  are  some  difference  and  some  additional  assumptions  for  MSDD. 

Consider  a  wireless  communications  system  operating  over  L  independent  reception 
branches,  the  sampled  output  of  the  matched  filter  corresponding  to  time  k  and  the  £-th 
branch  is 

rk,i  =  VPsetSk  +  zk,t ,  £  =  1, 2,  •  •  •  ,  L,  (5.1) 

where  Ps  is  the  power  of  the  desired  signal,  c/  is  the  channel  gain  of  the  £-th.  branch,  sk  is 
the  transmitted  M-DPSK  symbol,  and  zk,t  is  Gaussian  correlated  noise. 

The  transmitted  signals  can  be  expressed  as  sk  =  e]6k,  0k  =  27r(4-l)/M,  ik  - 
1,2,  •  •  •  ,  M.  The  transmitted  symbols  are  differentially  encoded,  i.e.,  6k  =  0*-i  +  A0*,  where 
A 9k  is  the  phase  representing  the  transmitted  information  at  time  k. 

The  signal  model  in  vector  notation  is 

rk  =  VPacsk  +  zk,  (5-2) 

where  vk  =  [r^  i,  *  *  *  ,  r  and  zk  are  vectors  defined  similar  to  Tk. 

The  channel  gains  c/’s  are  assumed  to  be  independent  and  identically  distributed  (i.i.d.), 
zero-mean,  circularly  symmetric,  complex  Gaussian  random  variables  (Rayleigh  fading),  with 
variance  Qt/2  per  dimension.  The  correlated  noise  term  zfc  is  the  aggregate  of  an  interference 
source  and  AWGN  and  it  is  assumed  to  be  complex- valued,  zero-mean,  circularly  symmetric, 
and  governed  by  a  Gaussian  distribution  with  covariance  matrix  ILZ  =  E  [z*zjf]  •  For  a  single 
interference  source  and  AWGN ,  the  covariance  matrix  can  be  expressed  as 

Rz  =  E  [zfczf  ]  =  P/C/cf  +  diag  (o\,  o\,  •  •  •  ,  <r|) ,  (5.3) 
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where  Pi  is  the  interference  power,  c j  is  the  interference  channel  vector,  and  (cr? ,  erf,  •  •  •  ,  a\) 
is  the  power  profile  of  the  AWGN. 

Consider  a  sequence  of  K  symbols  running  from  time  k  —  (K  —  1)  to  k.  Assume  the 
channel  is  static  over  the  duration  of  this  sequence.  Using  vector  notation, 

Lk  =  >/Ps Ms*,  +  zjt,  (5.4) 

where  r*  —  •  •  •  ,  r*]  ,  s*  —  •  •  •  , S/t]  and  is  a  vector  defined  similar 

as  r*.,  and  H  =  1K  <g>  c  is  the  channel  matrix  for  the  signal  of  interest,  where  ®  denotes  the 
Kronecker  product,  and  I K  is  the  identity  matrix  of  rank  K. 


5.3  Decision  Statistic 

We  formulate  the  decision  statistic  for  a  symbol  sequence  s*  =  •  •  •  ,sk]T  based  on 

an  observation  interval  consisting  of  length  K  as  embodied  in  the  vector  r*.  Assume  covari¬ 
ance  matrix  of  the  interference  plus  noise  R2  is  known,  the  maximum  likelihood  detector  for 
the  sequence  s*  is  given  by 


%  =  arg  max  p(rk\sk, R*),  (5.5) 

Sit 

where  p(rJSfc,R2)  is  the  likelihood  of  the  observed  data  rfc  given  the  transmitted  symbol 
sequence  sk  and  the  covariance  matrix  Rz.  Under  the  Gaussian  assumption  for  the  aggregate 
of  interference  and  noise,  the  observation  r*  conditioned  on  the  transmitted  sequence  s*, 
the  covariance  matrix  Rz  and  channel  c  has  a  multivariate  Gaussian  distribution.  The 
conditional  probability  p(rJt|sfc,Rz,c)  can  then  be  expressed  as 

p(r*|sfc,R*,c)  =  ir~KL\Rz\~K exp  j-  5^(r*-i  -  \Z^csfc_i)ffRj1(rfc_j  -  v^cs*_i)  j 

(5.6) 

Diagonalize  the  interference  plus  noise  covariance  matrix  R2  as  R2  =  U2A2Uf ,  where 
A2  =  diag  (Ai,  •  •  •  ,  AL) ,  Ai,«  •  •  ,A L  are  the  eigenvalues  of  R2,  and  U2  is  a  unitary  matrix 
whose  columns  are  the  eigenvectors  of  R2.  It  follows  that  (5.6)  can  be  written  as 

p(rfc|sit,R2,g)  =  exp  j-  sk-i)H h~l{xk-i  -  v/^gSfc-i)| , 

(5.7) 

where 

xfc_i  =  Ufr*_i  (5.8) 

g  =  Ufc.  (5.9) 

x*_j  is  the  whitened  received  signal  vector  and  g  is  the  modified  channel  vector.  Note  that 
since  U2  is  unitary,  the  modified  channel  vector  g  has  the  same  distribution  as  the  original 
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channel  vector  c.  Let  the  components  of  the  modified  channel  vector  g  be  expressed  as 
gt  =  a(ej<t>t ,  l  —  1, . . .  ,  L.  Likewise,  let  the  £th  component  of  xfc_j  be  xk-i,t.  Expanding  the 
exponent  in  (5.7)  and  grouping  terms  that  do  not  depend  on  g  or  sk-i,  we  obtain 

p(rJsfc,Rz,g)  =  v~kl\Rz\~k exp  {—Co} 

fl  exp  {-KPXtXo?t  +  2y/PsXjl  |  yt  (sfc)|  cos  {fa  -  0t  (s*))}  , 

1=  1 

(5.10) 


where 


K- 1  L 


c'«=EEv  i^-.  ,<r 

t=0  /=1 


(5.11) 


and 


(5.12) 


Note  that  yt  (s*)  is  a  function  of  both  the  transmitted  sequence  sk  and  the  observed  sequence 
Lk- 

Recalling  that  the  components  of  the  modified  channel  vector  g  have  the  same  distribution 
as  the  components  of  the  channel  vector  c,  it  follows  that  oti  is  Rayleigh  with  E  [o:|]  =  fl/  and 
fa  is  uniformly  distributed  in  the  interval  [0,  2i r).  To  average  the  conditional  distribution 
p(ij.|sfc,R2)g)  over  the  modified  channel  g,  we  need  to  evaluate  the  integral 


p(rifc|sJb,Rz)=  J p(r*|sfc,Rz,g)pg  (g)  dg  (5.13) 

where  pg  (g)  is  the  probability  density  function  of  g.  Assume  the  channels  are  independent 
to  each  other,  then 


L 

Ps  (g) = n  p*t  (&)  (5i4) 

t=i 

where  pae(at)  and  p<pl  (fa)  are  the  probability  density  functions  of  at  and  fa  respectively.  For 
Rayleigh  fading  channels, 


=  |^exp(-?0  0<a,<oo  (5.15) 

^  0  <  *  <  2tt.  (5.16) 
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Substitute  (5.10),  (5.14),  (5.15)  and  (5.16)  into  (5.13),  and  separate  the  integrations, 


I0{2y/FsXt  1 1  yt  (sk)\a()aedae, 


(5.17) 


p(rfc|sfc,Rz)  =  7T  km|Rz|  ^expl-Co}  exp  j-  ^PSA/  +  aedat 

exp  ^2\/PaXj1 1 yt  (sfc)|  at  cos  (<f>t  -  ^)}d^| . 

After  the  averaging  over  the  uniform  distribution  of  <f>t  is  carried  out, 
p(rk\sk,TLz)  =  Tr^lR^l^expf-Co} 


(5.18) 

(5.19) 


where  I0(x)  is  the  zeroth  order  modified  Bessel  function  of  the  first  kind.  Using  the  integra¬ 
tion  expression  of  Bessel  function  Iq  ( x )  in  Appendix  F,  we  obtain, 


p(r*:|sfc,Rz)  = 


7r  kl|Rz|  k exp  {-C0} 


KP&t  +  Xt) 


^<l^(sfc)l2  1 

Xe  ( KPsQt  +  A/)  J 


(5.20) 


In  (5.20),  only  the  argument  of  the  exponential  function  is  dependent  on  the  transmitted 
sequence  sk  since  only  the  terms  yt  (sjt)’s  are  functions  of  sk.  Due  to  the  monotonicity  of  the 
exponential  function,  maximizing  p(r*|s*)  with  respect  to  Sjt  is  equivalent  to  maximizing  the 
following  decision  statistic: 


t  \  — 

I,  Si  ”  ti  ^(kp,q,+x,y 


(5.21) 


From  (5.5),  the  corresponding  MSDD  decision  rule  is 

sfc  —  arg  max  r](sk).  (5.22) 

S  k 

The  detector  searches  through  sequences  s*  and  chooses  the  sequence  that  has  the  largest 
decision  metric  r/( sk).  A  diagram  of  the  MSDD  receiver  is  shown  in  Fig.  5.1. 


From  the  previous  relation,  it  follows  that  the  optimum  multiple  symbol  differential 
detector  for  multiple  channel  branches  and  in  the  presence  of  interference,  is  a  weighted  sum 
of  correlations  of  whitened  observations  and  hypothesis  symbols.  Note  that  this  decision 
statistic  does  not  require  knowledge  of  the  signal  channel  vector. 
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Si 


Figure  5.1:  The  diagram  of  multiple  symbol  differential  detector.  For  M-DPSK  with  M 
symbols,  the  number  of  symbol  sequences  that  need  to  be  tried  is  M\  =  MK 


The  decision  statistic  in  (5.21)  provides  multiple  symbol  differential  detection  for  a  M- 
DPSK  sequence  transmitted  over  multiple,  independent  fading  channels  in  the  presence  of 
correlated  Gaussian  noise. 

The  decision  statistic  is  ambiguous  with  respect  to  an  arbitrary  phase  O'.  Indeed,  let 
s'k  =  ej6'sk,  then 

|y<(s'fc)|  =  |y*(ei<’'sit)| 

=  ^2  Xk-i,i  (e3d'sk- i) 

=  \yi(sk)\  •  (5-23) 

Differential  encoding  at  the  transmitter  is  required  to  resolve  this  ambiguity. 


5.3.1  Iterative  Decision  Feedback  Algorithm 

The  complexity  of  MSDD  for  M-DPSK  with  a  K  symbols  observation  interval  increases 
with  MK~l.  For  large  K,  this  makes  simulations  impractical.  To  overcome  this  difficulty,  a 
practical  sub-optimal  algorithm  that  uses  decisions  feedback  was  implemented.  The  basic 
idea  of  the  algorithm  is  to  make  symbol  by  symbol  decisions  rather  than  testing  the  full 
sequence  of  symbols  simultaneously.  The  algorithm  proceeds  from  symbol  to  symbol  along 
the  sequence  of  K  symbols;  at  symbol  i  it  maximizes  a  decision  statistic  assuming  that  the 
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Figure  5.2:  Comparison  of  optimum  algorithm  and  iterative  decision  feedback  algorithm  for 
L  =  4  branches,  DPSK  modulation,  SIR  =  —6  dB. 


other  ( K  —  1)  symbols  have  been  detected  and  are  known.  Several  iterations  can  be  carried 
out  to  improve  performance.  The  algorithm  was  implemented  as  the  following  procedure: 

1.  Initialization: 


(a)  Initialize  iteration  index,  m  =  0. 

(b)  Initializes^.  = 

(c)  Initialize  time  index  i  =  K  —  2. 


2.  Increase  iteration  index  m  +  1  — ►  m. 

3.  For  i  =  (K  —  2)  to  0, 


Evaluate 
End  loop  i. 


arg  max  rj 

a(m) 

Sk-i 


([ 


-l  (m) 

Sk-(K-2)i ' 


(m)  (m— 1) 


(m— 1) 


4.  If  m  is  not  equal  to  the  required  iteration  number  (which  is  determined  empirically), 
go  back  to  step  2. 

5.  Differentially  decode  sj^  to  get  the  final  output. 
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To  demonstrate  the  performance  of  this  sub-optimal  decisions  feedback  algorithm,  Fig. 
5.2  compares  the  sub-optimal  and  optimal  (based  on  (5.22))  algorithms.  The  comparison 
is  for  the  case  of  L  =  4  diversity  branches,  DPSK  modulation,  and  SIR  =  —6  dB.  For  an 
observation  interval  of  K  —  12  symbols,  with  just  2  iterations,  the  performance  of  the  sub- 
optimal  decisions  feedback  algorithm  is  within  just  0.2  dB  of  that  of  the  optimum  algorithm. 
The  advantage  of  the  sub-optimal  decisions  feedback  algorithm  is,  of  course,  that  it  takes 
much  less  time  to  run  than  the  optimum  algorithm.  From  the  figure,  it  can  also  be  observed 
that  iterations  are  beneficial  to  the  performance  of  decisions  feedback.  The  second  iteration 
provides  about  0.5  dB  gain  relative  to  that  of  without  iteration  (iteration  1).  Additional 
iterations  do  not  seem  to  improve  the  performance. 


5.3.2  Special  Case 

Some  special  cases  provide  insights  into  the  operation  of  MSDD.  For  a  channel  with  a  flat 
gain  profile  f 1/  =  1,  and  a  flat  AWGN  profile  <r2  =  <r2  for  i  =  1,2, •••  ,L,  (5.21)  can  be 
expressed  as 


L 


v(*k)  =  53 

t-l 


|y<(s«:)l2 

A,  (KP.  +  XiY 


(5.24) 


We  further  specialize  (5.24)  to  the  following  special  cases. 


No  Interference 

For  this  case,  Pj  =  0,  the  noise  covariance  matrix  Rz  =  a1! L,  eigenvalues  A e  =  a2,  £ 


1, 2,  •  •  •  ,  L,  and  Uz  =  I L.  Then  (5.12)  simplifies  to 

K- 1 

Vi  (s*)  =  J3  rk-i,t4-i- 
»= 0 

The  decision  statistic  in  (5.24)  becomes 

1  r  2 

=  cr2  {KPS  +  <72) 


(5.25) 


(5.26) 


Since  the  term  outside  the  sum  is  independent  of  s*,  the  above  decision  statistic  is  equivalent 
to 

|2 


*?(s  k)  =  '52\yt(.sk)\2  =  ^2 


1=  1 


fol 


K- 1 


53  r^sk~i 


i—0 


(5.27) 


This  decision  statistic  is  the  same  as  that  in  [31,  equation  (8)].  Indeed  (5.21)  is  the  general¬ 
ization  of  [31,  equation  (8)]  to  MSDD  in  the  presence  of  interference. 
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Interference  Noise 


For  a  uniform  AWGN  power  profile,  the  eigenvalues  of  the  interference  plus  noise  covariance 
matrix  (5.3)  are  Ai  =  Pj  £Li  |c/,*|2+a2,  A2  =  . . .  =  AL  =  a1-  For  a  high  interference  to  noise 
ratio,  Ai  »  Xt,  l  ^  1.  It  follows  that  the  decision  statistic  in  (5.24)  can  be  approximated  by 
the  expression 


/  \  |j/l(sfc)| 

V S  (■ KPs  +  v 2Y 


(5.28) 


The  interpretation  of  this  result  is  that  for  a  strong  interference  source,  the  decision 
statistic  is  similar  to  that  of  MSDD  without  interference  and  one  fewer  degree  of  freedom. 
This  result  will  be  further  demonstrated  in  the  ensuing  error  probability  analysis. 


5.4  Error  Probability  Analysis 

An  exact  expression  for  the  BEP  for  differential  detection  can  be  obtained  only  for  DPSK 
modulation  and  K  =  2  symbols.  The  exact  error  analysis  is  intractable  for  the  general 
case  of  MSDD  with  M-DPSK  modulation  over  diversity  channels  and  in  the  presence  of 
interference.  The  alternative  approach  is  to  obtain  an  analytical  approximate  upper  bound. 
In  this  section,  we  first  derive  an  exact  expression  for  the  PEP.  Then,  using  this  expression, 
we  derive  the  union  bound  of  the  BEP.  From  the  union  bound  an  approximate  upper  bound 
is  derived.  The  approximate  upper  bound  consists  of  relatively  simple  algebraic  expressions. 
Even  simpler  expressions  are  obtained  for  the  asymptotically  large  SNR  and  small  signal  to 
interference  ratio  (SIR).  In  the  numerical  results  section,  it  is  shown  that  the  approximate 
upper  bound  is  very  close  to  the  BEP  obtained  by  simulation. 

5.4.1  PEP  Analysis 

In  the  derivation  of  the  PEP,  we  assume  a  uniform  flat  power  profile  for  the  desired  signal 
channel,  Qt  =  1,  and  a  flat  AWGN  profile  with  a\  =  <r2  for  l  =  1,2,  •  •  •  ,  L.  The  PEP  is 
developed  for  correlated  noise  characterized  by  the  covariance  matrix  in  (5.3). 

In  general,  the  interference  source  is  subject  to  effects  of  the  fading  channel  (similar 
to  the  desired  source).  It  follows  that  analysis  using  the  covariance  matrix  R2  in  (5.3)  is 
conditional  on  the  interference  random  channel  c /.  Results  obtained  from  such  analysis  need 
to  be  averaged  over  the  distribution  of  Cj.  Fortunately,  this  complication  can  be  avoided 
recognizing  that  when  the  detector  acts  to  suppress  the  interference,  there  is  only  a  small 
penalty  in  using  in  the  analysis  the  average  value  of  the  interference  power  PjE  jcf  Cj]  in 
lieu  of  the  instantaneous  power  P/cf  Cj  (see  [15]).  Assuming  that  the  interference  channel 
c /  is  complex- valued,  zero  mean  and  with  variance  fI/^/2  =  1/2  per  dimension,  it  follows 
that  the  average  eigenvalues  of  Rz  are 


Ai  =  PiE  [cf  C/]  +  <r2 
=  LP!+a2 


(5.29) 
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and  A*  =  <x2  for  l  =  2, 3,  •  •  •  ,  L. 

Let  s k  and  s'k  denote  two  sequences  each  containing  K  M-DPSK  symbols.  The  PEP  that 
Sfc  is  transmitted  but  s*  is  detected  (sfc  7^  Sosfc>  where  So  is  an  arbitrary  M-PSK  symbol)  is 
denoted  as  P  (s*  — >  sj.) .  An  error  event  occurs  when  ^(s*)  <  r?(s'fc).  Define  random  variable 

D, 


D  =  f?(s k)  -  v{s'k). 


(5.30) 


Note  that  D  is  random  due  to  both  the  random  noise  and  the  random  channels.  We  seek  to 
evaluate  the  probability  that  D  <  0. 

Using  steps  similar  to  [7,  Appendix  B],  it  can  be  shown  that 


P  (Sfc  ->  s'*) 


_l  r+i* 

2 Tfj  J — oo+j« 


=  -  ^  Res 


Im(ai/)>0 


MM. 


U) 


\Ul 


(5.31) 


where  e  is  a  small  positive  number;  is  the  characteristic  function  of  D\  Res  u)t 

denotes  the  residue  of  at  pole  uf,  the  summation  is  taken  over  the  poles  in  the  upper 

half  of  the  complex  plane. 

In  Appendix  G,  the  following  expression  is  derived  for  the  characteristic  function  of  the 
random  variable  D : 


$r>(jw)  ^  ^  _  j^w)  ^  _  j^3 W)L  1  (1  _  jniU>)L  1  ’ 


(5.32) 


where 


m  = 


\b\  (C P.  ±  VW2  +  4  (KPS  + Ax)  CAx) 
\b\  [CPs  ±  y/?P*  +  4  (A'P,  +  A2)  CA2) 


*  =  1,2 

i  =  3,4 


(5.33) 


and 


/  r_~ 

(5.34) 

V  Ai  ( KPa  +  Ai) 

/  i”- 

(5.35) 

V  A2  (KPS  +  a2) 

K2  -  \v(sk,s'k)f . 

(5.36) 

In  (5.33),  the  plus  sign  is  taken  for  i  =  1  and  3,  while  the  minus  sign  is  taken  for  i  =  2  and 
4.  In  (5.36),  v(sk,s'k)  =  skffsk  is  the  correlation  coefficient  between  the  transmitted  sequence 
sk  and  the  detected  sequence  s'k.  Note  that  0  <  v(sk,  s'k)  <  K,  with  equality  on  the  right 
hand  side  when  s*  =  s*. 
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Substituting  (5.32)  into  (5.31),  the  PEP  is  obtained  as 
P  (sfc  -4  s'fc) 


=  —  ^  Res 

lm(ui,)>0 


W  (1  -  JMi^)  (1  -  J>2w)  (1  -  j/i3w)  (1  -  jn 4w) 


(5.37) 


Since  fj.i,  fj,3  >  0  and  //2,  // 4  <  0,  only  the  poles  -j//J2  and  —  j / /x4  are  in  the  upper  half  of 
the  complex  plane.  Eq.  (5.37)  becomes 


P  (sfc  s'k ) 

1  1 


=  —Res 
—Res 


1 


Lw  (!  “  iMi")  (!  “  Wiu)  (1  -  jfi3u))L  1  (1  -  jiMu) 
111  1  1 


1  .  1 
L— 1  ’  3 


[uj  (1  -  juxu)  (1  -  (1  -  j^u)L  1  (1  -  juiuj) 


M2J 

.  1 


(5.38) 


After  some  cumbersome,  but  straightforward  manipulations  (which  are  similar  to  the 
manipulations  shown  in  Appendix  B),  it  can  be  shown  that  the  PEP  is 


P  (s/t  -t  s'*) 


(M2) 


21,-1 


(-i)Vr1  i 


L- 2 


(M i  -  M2)  (m 3  -  M 2)  (M4  -  M2) 

r  Ml 


L— 1 


Ms  1 


(£-2)'S 


(L-2  +  fc)! 

k\ 


(M 1  “  M2)  (^1  -  /x4) 


Mi _ + _ M 2 


M^1-* 


L-l-k 


(Mi  ~  M2)  (M2  -  M4) 


L-l-fc 


^-i+* 


(«>  -  mi) 


t-i+fc • 


(5.39) 


The  former  expression  is  the  exact  PEP  of  MSDD  with  diversity  branches  and  a  rank  one 
interference  source.  The  PEP  is  conditioned  on  the  transmitted  sequence  s*  and  is  a  function 
of  the  detected  sequence  s'fc.  Note  that  the  expression  in  (5.39)  already  incorporates  statistical 
information  on  the  channel  and  interference.  This  form  of  the  PEP  is  quite  complicated  and 
does  not  afford  much  insight.  It  is  of  interest  to  obtain  simpler  expressions  for  special  cases. 
In  the  ensuing  analysis,  the  symbol  signal  to  noise  ratio  SNR  is  7  =  Pa/ a2,  and  signal  to 
interference  ratio  is  SIR  =  Ps/Pi- 


No  Interference 

For  this  case  Pj  =  0,  Ai  =  A2  =  a2.  From  (5.33)  we  have: 


Mi  = 


_! _ lL_ 

2  P,  (*7+1) 

1 

2 Ps  (K7+ 1) 


/C!  +  4 

K  +  { 

/c2  +  4 

K  +  i 

i  =  1,3 


i  =  2,4 


(5.40) 
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(5.41) 


Since  Hi  =  /x3,  and  /*2  =  /i4,  it  follows  that  the  PEP  in  (5.38)  can  be  rewritten  as 


P  (Sjfc  -4  sjfe)  =  -Res 


'1  1 _ 1  .  1 ' 

v  (1  -jniu)L  (1  -JH2u)L'  3 M2_ 


It  is  not  too  difficult  to  show  that  this  results  in  the  following  expression  for  the  PEP: 


P  (s*  -4  s'k)  = 


(Mi  ~  M2 )l(L  -  1)! 


yv  (L—1  +  l)!  Hi 


L  £—0 


t>-  (mi  -  Ml)' 


(5.42) 


It  could  be  shown  numerically  that  (5.42)  is  equivalent  to  the  PEP  developed  in  [31].  How¬ 
ever,  (5.42)  has  the  advantage  that  it  provides  the  PEP  in  closed  form  without  the  need  of 
integration. 

The  case  of  no  interference  can  be  further  simplified  for  large  signal  to  noise  ratio  7  »  1. 
In  this  case,  (5.40)  simplifies  to 


„(7  C/KPs  <  =  1,3 
-1/P,  *  =  2,4  • 


Mi 

Substituting  these  results  in  (5.42)  and  noticing  that  hi  ^  M2>  we  have 

Pfs  .sn  ~  (1/P,)L  yifr-i  +  Qi 

(*  k)  ('rC/KPs)L(L- 1)!^  <! 

'2L  -  1\  1 

L  J(h)L' 

This  expression  clearly  exhibits  the  L-order  diversity  of  the  system. 


c 


(5.43) 


(5.44) 


SIR  <  1,  SNR  >  1 

By  assumption  Pj  »  Pa  a2,  therefore  Ai  Ps  a2.  From  (5.34),  we  have  61  ~  1/Ax. 
After  some  simple  manipulations,  it  follows  from  (5.33), 

Mi, 2  *  (5.45) 

To  evaluate  Hz  and  //4,  we  approximate  62  ~  \j\jKPao'1  and  substitute  in  (5.33)  to  obtain 
Hz  ~  C/  (A" cr2)  and  Hi  ~  —  1/P,.  Substituting  these  approximate  values  into  (5.39)  and 
keeping  only  the  dominant  term,  after  some  manipulations,  we  have 

(546) 

Comparing  (5.46)  with  (5.44),  we  can  see  the  PEP  for  systems  with  diversity  L  and  a  large 
interference  is  equal  to  the  PEP  for  systems  with  diversity  ( L  —  1)  and  without  interference. 
This  result  is  well  known  for  interference  suppression  using  OC.  This  analysis  proves  that 
the  loss  of  degree  of  freedom  due  to  interference  suppression  carries  over  to  MSDD  over  a 
diversity  Rayleigh  channel. 
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5.4.2  BEP  Approximate  Upper  Bound 

The  sequence  sfc  of  M-DPSK  symbols  corresponds  to  (K  -  1)  log2  M  information  bits  (with 
differential  encoding,  the  first  symbol  is  known).  Let  ufc  be  the  sequence  of  (K  -  1)  log2  M 
information  bits  encoded  as  Sjt,  and  let  be  the  sequence  of  information  bits  which  results 
from  the  detection  of  s'k.  The  pairwise  BEP  associated  with  transmitting  a  sequence  ufc  and 
detecting  another  sequence  uj.  is  given  by 

Pb  (sfc  -4  s'k)  =  ( Sk  sfc)  ’  (5‘47) 

where  h(uk,  uj.)  denotes  the  Hamming  distance  between  U*  and  n'k. 

The  BEP  that  sk  is  transmitted,  but  an  error  sequence  (any  error  sequence)  is  detected, 
is  upper  bounded  by  the  union  of  all  pairwise  bit  error  events.  Since  s*  can  be  any  input 
sequence  (e.g.,  the  null  sequence  sk  =  [1, 0, 0,  •  •  •  ,  O]7'),  we  drop  the  dependency  on  s*  from 
the  notation.  The  union  bound  on  the  BEP  can  then  be  written  as 

Pb  < 


where  the  summation  is  taken  over  all  the  sequences  u*.’s  which  are  different  from  the  trans¬ 
mitted  sequence  of  information  bits  u*. 

Direct  application  of  (5.48)  does  not  shed  light  on  the  mechanisms  affecting  MSDD 
performance.  A  clearer  picture  is  obtained  by  developing  an  approximation  to  the  union 
bound.  Note  that  the  union  bound  in  (5.48)  is  a  function  of  the  PEP’s,  which  in  turn  are 
determined  by  H2,  H3  and  A4  (see  (5.39)).  From  (5.33),  /xi,  /X2,  /J3  and  Hi  are  functions  of 
the  quantity  |v(sfc,  s*)|2  through  the  relation  (  =  K2  -  |v(sfc,  s'fc)|2 .  In  [24],  it  is  shown  that 
on  the  AWGN  channel,  for  large  SNR,  the  dominant  terms  in  the  BEP  occur  for  sequences 
for  which  the  quantity  |t;(sfc,  s*)|2  is  maximum.  Carrying  over  the  same  approach  to  the 
fading  channel,  keeping  only  the  dominant  terms  and  noticing  that  P  (sk  — >  s*)  is  constant 
if  \v  (sfc,  s*)|  is  constant,  we  obtain  the  following  approximation  to  the  union  bound 


Y  Pb(Sk  S*) 


(jr-i)iS£5 


(5.48) 


A  = 


h{  ufc,  u'fc) 


(K  —  1)  log2  M 

I 

|v(sfc,s'fc)|  =  |w(sfc,s'fc) ! 

[p(s*  -t  «yi|.(My,)H.(«y,)|  J  ■ 

The  maximum  value  of  |u(sfc,s'fc)|  was  shown  in  [24,  Eq.  (38)]  to  be 


l«  K,  =  ^(A:- lf  +  2(/f-l)( 1  -  2  sin2  ^  )  +  1. 


(5.49) 


(5.50) 
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Also  from  [24,  Appendix  B],  for  sequences  such  that  |v(sj|.,sjj. 
mutated  Hamming  distances  are 

£  ■*(«*."»)  =  {  2(jr-i), 

“ik/Ufc 

for  binary  modulation,  M  =  2  and 

X)  MU*X)  =  {  4(^1 1); 

u'fc#U)fc 

|V(Sfc,s'fc)|  =  |v(sfc>s'fc)|max 

for  multilevel  modulation,  M  ^  4. 

Strictly  speaking,  (5.49)  is  not  an  upper  bound  of  the  BEP.  Numerical  results,  however, 
show  that  it  is  very  close  to  or  larger  than  the  BEP  obtained  by  simulation.  Therefore  we 
will  use  (5.49)  to  study  the  performance  of  MSDD  in  the  presence  of  interference. 

Next,  we  evaluate  the  approximate  upper  bound  for  differential  binary  PSK  (DPSK)  and 
M-DPSK  (M  >  4)  modulations. 

DPSK  (M  =  2) 

For  this  case,  from  (5.50)  we  have 

=  (553> 

For  conventional  differential  detection,  the  observation  interval  is  K  =  2  symbols,  |u  (s*,  s*)|max 
0.  In  this  case,  there  is  only  one  error  sequence,  therefore  the  PEP  is  also  the  BEP, 

Pt=P(ss-tsi)|w„yi)H.  (5.54) 

Substituting  P(s*  -4  s'k)  from  (5.42)  into  (5.54),  we  obtain  the  exact  BEP  for  DPSK 
over  L  diversity  fading  channels  without  interference.  For  high  SNR  1,  using  (5.44),  we 
get 

A.DPSK  (“j.  ‘)p^E-  <6-65) 

This  expression  is  the  same  as  the  one  in  [7,  Eq.  (14-4-28)],  and  it  demonstrates  that  familiar 
expressions  for  differential  detection  can  be  obtained  as  a  special  case  of  the  general  case 
treated  in  this  report. 

For  a  longer  observation  interval  K  >  2,  substitute  (5.53)  and  (5.51)  into  (5.49)  to  obtain 

^dpsk  =  2  P  (sj;  — >■  s/fe)||u(SfciS»j|=fir_2  •  (5.56) 

This  expression  is  for  the  approximate  BEP  upper  bound  for  DPSK  over  slow-fading  Rayleigh 
diversity  channels  with  interference.  Next,  we  compute  some  special  cases  for  K  >  2  and 
SNR  7  >•  1,  which  result  in  simplified  expressions. 


)l  =  |»(Sfc,  sjfc)lmax’  the  aCCU- 


K  =  2 
K  >  2 


(5.51) 


K  =  2 
K  >  2 


(5.52) 
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No  Interference,  SNR  »  1  Using  (5.44)  in  (5.56),  the  approximate  upper  bound  is 


-d-DPSK 


|»(gfe,s)fe)|=K--2 


2^-^ _ l _ . 

V  L  )  4 L^L  (l  _  l.)L 


(5.57) 


SIR<  1,  SNR  >  1  Substituting  (5.46)  in  (5.56)  we  obtain  a  result  similar  to  (5.57), 
except  L  is  substituted  with  (L  —  1) , 


-Adpsk 


/2  (L  —  1)  —  1\  1 

V  L-l 


(5.58) 


M-DPSK  (M  ^  4) 

For  M-DPSK,  substitute  (5.52)  and  (5.50)  into  (5.49),  we  obtain  the  following  approximate 
upper  bound: 


^4m-DPSK  =  -+  S*)l|»(sfclS'fc)|=2|cos^| 


(5.59) 


for  K  =  2  symbols  and 

Am‘dpsk  =  P  (Sfc  ^  ^ 1  ><)  |=v'(K-x)i+2(Jir-i)(i-2si^  £)+i 

for  observation  intervals  of  size  K  >  2. 

Simplified  expressions  for  special  cases  are  computed  below. 


(5.60) 


No  Interference,  SNR  >  1  Using  (5.44),  we  obtain  from  (5.59)  and  (5.60)  the  approxi¬ 
mate  upper  bound 


2  f2L— 1\  1 

^M-DPSK  -  \^M  \  L  )  2  V  sin2L  (t r/M) 


(5.61) 


for  K  =  2  and 

,  4  (2L-l\ _ 1 _ 

m  dpsk  i0g2  ^y  X,  J  4  LjL  gin21-  (71-/M)  (l  — 


(5.62) 


for  K  >  2. 
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SIR  <C  1,  SNR  1  Substituting  (5.46)  in  (5.59)  and  (5.60)  we  obtain  a  result  similar  to 
(5.61)  and  (5.62),  except  that  L  is  substituted  with  (L  -  1) , 


.  2  /2(L-l)-l\ _ 1 _ 

m-dpsk  L  —  l  )  2i~17L-1  sin2^-1)  ( tt/M ) 


(5.63) 


for  K  =  2  and 

,  __J_/2(L-1)-1\ _ I _ 

m-dpsk  i0g2  M  \  L  —  l  )  \L-iryL-i  sin2(L_1)  (tt/M)  (l  —  j^)L 


(5.64) 


for  K  >  2. 


5.5  Comparison  with  OC 

In  this  section,  we  compare  the  BEP  of  MSDD  with  that  of  OC  analytically.  Since  the  expres¬ 
sions  of  the  BEP  for  general  cases  are  very  complex,  we  can  only  compare  the  performance 
analytically  for  the  cases  of  small  SIR  (relative  large  interference)  and  no  interference.  Since 
both  cases  yield  similar  results,  we  only  do  that  for  the  case  of  small  SIR.  In  next  section, 
we  will  compare  the  performance  of  MSDD  and  OC  numerically. 

From  (4.25)  and  (4.26)  in  Chapter  4,  for  SNR  7  >  1  and  SIR  <  1,  the  BEP  for  OC  is 
approximated  by  the  expressions 

p  ~2l'2(i_1)~1i  — _  (5.65) 

-Pb,BPSK  ~  2  ^  (t  _  !)  _  l  )  4L-171-1  1  > 


and 


4  (2  (L-  1)  -  1\ _ 1 _ 

•Pb,M-PSK  ~  \0g2M  \  (L  —  1)  —  1  /  i^-iryL-i  sin2(L_1l  (tt/M) 


(5.66) 


The  above  two  equations  are  for  OC  with  M-PSK  modulation.  Since  MSDD  uses  M-DPSK 
modulation.  We  want  to  compare  them  on  the  same  basis. 

For  OC,  the  exact  expression  of  the  BEP  for  M-DPSK  is  very  difficult  to  obtain.  But 
judging  from  Eq.  (4.200)  in  [37],  the  BEP  for  M-DPSK  is  about  twice  that  of  M-PSK 
except  for  very  small  SNR.  That  can  be  demonstrated  by  simulation.  Therefore,  the  BEP 
for  M-DPSK  is 


Pb, M-DPSK  ~  2 Lb, M-PSK- 


(5.67) 


Substitute  (5.65)  and  (5.66)  in  (5.67),  we  obtain  the  BEP  for  OC  using  M-DPSK  as 


(2  (L-  1)  -  1\  1 

Fb,DPSK«2^(L_1)_1 


(5.68) 


and 


4  (2  (L  —  l)  —  1^ _ 1 _ 

Lb, m-dpsk  ~  i0g2ji/  y  (L  —  1)  —  1  )  4^-i7^-1  sin2(L-1^  (tt/M) 


(5.69) 
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The  ratio  of  the  BEP  of  OC  and  the  approximate  upper  bound  of  MSDD  (for  K  >  2)  is 
given  by  the  ratios  of  (5.68)  to  (5.58)  and  (5.69)  to  (5.64),  respectively 

fi>,DPSK  _  -Pb,M-DPSK 
Adpsk  -^M-DPSK 

The  ratio  approaches  1  as  K  — t  oo.  We  conclude  that  in  the  presence  of  interference,  when 
the  observation  interval  of  MSDD  increases  to  infinity,  i.e.,  K  oo,  the  performance  of 
MSDD  approaches  that  of  OC  with  differential  encoding. 

According  to  (5.67),  the  BEP  of  OC  with  differential  encoding  is  about  twice  that  of  OC 
without  differential  encoding.  Therefore  for  MSDD  of  large  K,  the  BEP  is  only  about  twice 
that  of  OC  without  differential  encoding. 


-a 


L—l 


(5.70) 


5.6  Numerical  Results 

Numerical  results  presented  in  this  section  include  Monte  Carlo  simulation  results  and  anal¬ 
ysis  results.  In  all  cases,  the  channel  branches  and  noise  power  profiles  are  assumed  to  be 
uniform,  i.e.,  Q/  =  1  and  o\  =  <j2  for  i  =  1, 2,  —  ,  L.  The  bit  SNR  is  (Ps/log2(M))  /c2.  For 
comparison  purposes,  we  also  provide  BEP  curves  for  OC  with  differential  encoding.  All  the 
figures  are  for  L  =  4  diversity  branches. 

Fig.  5.3  shows  the  BEP  versus  SNR  for  DPSK  at  SIR  =  -6  dB.  Curves  labeled  ‘Sim¬ 
ulation’  represent  simulation  results,  while  curves  labeled  ‘Analysis’  show  analytical  results 
as  yielded  by  the  approximate  upper  bounds  (5.54)  (for  K  —  2)  and  (5.56)  (for  K  >  2).  In 
all  cases,  PEP’s  were  exact  as  computed  by  (5.39).  The  interference  plus  noise  term  was 
generated  such  that  its  covariance  matrix  followed  (5.3).  The  OC  curve  was  generated  by 
simulation.  It  can  be  observed  that  analysis  results  are  very  close  to  simulation  results.  It 
is  also  observed  that  the  performance  of  MSDD  approaches  that  of  OC  with  differential  en¬ 
coding  as  K,  the  number  of  symbols  in  observation  interval  increases.  For  example,  at  BEP 
=  2  x  10-3,  when  K  =  2,  the  SNR  difference  between  MSDD  and  OC  is  about  2.2  dB.  When 
K  =  7,  the  difference  is  about  1.0  dB.  At  K  =  40,  the  difference  becomes  an  insignificant 
0.2  dB. 

Fig.  5.4  and  Fig.  5.5  are  for  DQPSK  and  8-DPSK  respectively.  The  curves  in  these 
figures  follow  the  same  trends  as  in  Fig.  5.3. 

The  results  shown  in  Fig.  5.6  to  5.8  are  all  analytical  results.  In  these  figures,  bit  error 
probabilities  are  represented  by  their  approximate  upper  bounds.  The  approximate  upper 
bound  is  computed  based  on  the  exact  PEP  expression  in  (5.39)  except  for  Fig.  5.8. 

Fig.  5.6  shows  the  BEP  of  MSDD  as  a  function  of  the  number  of  symbols  in  the  obser¬ 
vation  interval,  K.  It  is  evident  that  for  both  DPSK  (binary  modulation)  and  for  8-DPSK 
(M  =  8),  the  performance  of  MSDD  approaches  that  of  OC  as  the  observation  interval 
increases. 

Fig.  5.7  shows  the  BEP  versus  SIR,  for  bit  SNR  =  10  dB,  and  for  the  cases  of  K  =  2 
and  K  -  40  symbols.  It  is  observed  that  when  K  =  40,  MSDD  achieves  performance  close 
to  that  of  OC  with  differential  encoding  regardless  of  the  SIR. 

Fig.  5.8  is  intended  to  verify  the  asymptotic  large  SNR  approximation  to  the  PEP.  The 
signal  modulation  is  DQPSK.  Curves  labeled  ‘asymp’  represent  asymptotic  results  computed 
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by  applying  (5.63)  (for  K  =  2)  and  (5.64)  (for  K  =  40);  curves  labeled  ‘exact’  represent 
exact  results  from  (5.59)  (for  K  =  2)  and  (5.60)  (for  K  =  40).  It  is  observed  that  for  most 
SNR  of  interest  (SNR  >  10),  the  approximate  upper  bound  based  on  asymptotic  PEP  is 
very  close  to  the  approximate  upper  bound  based  on  the  exact  PEP. 
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Figure  5.3:  BEP  versus  SNR  for  L  —  4  branches,  DPSK  modulation,  SIR  =  —6  dB. 


Figure  5.4:  BEP  versus  SNR  for  L  =  4  branches,  DQPSK  modulation,  SIR  =  —6  dB. 
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SIR  (dB) 

Figure  5.7:  BEP  versus  SIR  for  L  —  4  branches,  DQPSK  modulation,  bit  SNR  =  10  dB. 


Bit  SNR  (dB) 


Figure  5.8:  Comparison  of  asymptotic  results  and  exact  results  for  L  =  4  branches,  DQPSK 
modulation,  SIR  =  —6  dB. 


Chapter  6 


Multipe-Symbol  Differential 
Detection  with  Known  Channel 
Information  of  Interference 

6.1  Introduction 

Until  now,  we  have  discussed  two  kinds  of  detectors  for  communication  systems  with  recep¬ 
tion  diversity  in  the  presence  of  white  Gaussian  noise  and  interference  source.  In  Chapter 
3,  we  discussed  detector  using  OC.  To  implement  OC,  the  channel  gain  of  the  desired  signal 
and  the  covariance  matrix  of  the  interference  plus  noise  must  be  available  to  the  receiver. 
In  Chapter  5,  we  discussed  MSDD  for  the  case  where  the  channel  gain  of  the  desired  signal 
was  assumed  to  be  unknown,  but  the  covariance  matrix  of  the  interference  plus  noise  was 
assumed  to  be  known.  Both  detectors  show  the  ability  to  suppress  interference. 

It  would  be  desirable  to  be  able  to  suppress  the  interference  without  requiring  any  infor¬ 
mation  of  the  interference.  Unfortunately  that  is  impossible.  In  Chapter  3,  the  interference 
plus  noise  is  modeled  as 

Ni 

Zfc  = 

*=i 

If  we  don’t  have  any  information  about  c*  (which  is  assumed  to  be  Gaussian  distributed) , 
the  interference  term  \fPi  YaI\  c iSi,k  would  be  the  same  as  Gaussian  noise  and  could  not  be 
distinguished  from  the  white  Gaussian  noise  n*.  Therefore  at  least  some  information  about 
the  interference  is  required. 

In  this  Chapter,  we  develop  detector  for  the  case  where  the  only  required  channel  in¬ 
formation  is  the  amplitude  of  the  channels  of  the  interference.  The  scenario  is  similar  to 
that  in  chapter  5.  But  in  addition  to  assuming  that  the  channel  gain  of  the  desired  signal 
is  unknown,  the  phase  of  the  channel  of  the  interference  is  assumed  to  be  unknown  as  well. 
The  channel  amplitude  of  the  interference  is  assumed  known.  Moreover,  the  interference  is 
assumed  to  have  the  same  MDPSK  modulation  as  the  desired  signal.  A  maximum  likelihood 
sequence  detector  (MLSD)  is  formulated  for  the  joint  detection  of  the  desired  signal  and  the 
interference.  Simulation  is  performed  for  DPSK  modulation.  Simulation  results  in  terms  of 


55 


BEP  versus  SNR  are  given  out  and  compared  with  the  results  obtained  by  other  detection 
schemes. 

It  is  shown  that  when  the  interference  level  is  high,  this  MSDD  technique  can  achieve 
better  performance  than  detectors  using  OC  (with  differential  encoding). 


6.2  System  Model 

The  system  model  used  in  this  chapter  is  similar  to  that  mentioned  in  Section  2.2,  except 
that  now  we  assume  there  is  only  one  interferer.  The  output  of  the  match  filter  is 

rk,t  =  y/PsCiSk  +  VPicIitSl>k  +  nkii,  l  =  1, 2,  •  •  •  ,  L.  (6.2) 

The  definitions  of  the  variables  are  in  Section  2.2.  Both  the  desired  signal  sk  and  the 
interference  source  sI>k  are  assumed  to  be  M-DPSK  symbols. 

The  signals  in  vector  notation  are 

rfc  —  \/PsCsk  +  JPlClS[,k  +  nk.  (6.3) 

We  assume  both  ct  and  cItt  are  zero-mean  complex  Gaussian  random  variables  (Rayleigh 
fading),  and  they  are  mutual  independent.  For  convenience,  we  define  q  =  aee^e,cn  = 
vectors  a  =  [c*i,  £*i,  •  •  •  ,  on]  •  And  similarly  as  a,  define  vector  <f>,  aj  and  <f>j.  In 
this  paper,  aT  is  assumed  to  be  known,  but  a,  <j)  and  <£/  are  assumed  to  be  unknown. 
Assume  the  covariance  matrix  of  nk  is 

Bm  =  E  [n*nf ]  =  diag  [aj,  of,  •  •  •  ,  of] ,  (6.4) 

where  aj  (£  =  1,  •  •  •  ,  L)  is  the  power  of  the  noise  on  the  Ath  branch. 

Consider  a  sequence  of  K  symbols  running  from  time  k  —  (K  -  I)  to  k.  Assume  the 
channels  are  static  within  the  duration  of  this  sequence.  Using  vector  notation, 

+  y/PiHjSiik  +  n*,  (6.5) 

where  r*  =  [r*_(^_1),rfc_(^_2),  •  •  •  ,r*]T.  The  channel  matrix  H  =  Ik<3>c,  where  <g>  denotes 
the  Kronecker  product,  and  lK  is  the  identity  matrix  of  rank  K.  H/  =  I^®C/,  s*  = 

[s*  (*-_!), ■■■  ,sk]  ,  s =  [s/,*— (*:— 1), •  •  •  , s/)fc]T ,  n*.  =  [nfc_(Ar_i), •  •  •  , nfe]r. 

Some  other  assumptions  of  the  signals  and  channels  will  be  given  out  in  the  derivation 
of  the  decision  statistic. 


6.3  Decision  Statistic 

In  this  section,  we  derive  the  decision  statistic  for  MSDD  with  known  channel  amplitude 
of  the  interference  aj.  MSDD  is  a  form  of  MLSD.  The  decision  is  made  after  K  symbols 
are  transmitted  and  received.  Conditioned  on  the  channels  c  and  c /,  the  coherent  decision 
criterion  for  sequence  detection  is  given  by 

(sk,sTtk)  =  arg  max  p(r*|sfc,  s r>k,  c,  c7).  (6.6) 

sfcjs/,fc  v  ' 
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Note  that  we  perform  detection  of  both  the  desired  and  the  interference  symbols.  The  pair 
(Sfc,  s itk)  that  maximizes  p(rfc|sfc,  s /,*,  c,  c j)  is  chosen  as  the  detected  symbols.  When  the  only 
known  channel  information  is  a/,  the  decision  rule  for  MLSD  would  be 


(sit,  s/)fc)  =  arg  max  p(r*|sfc, S/,fc,  <*/)•  (6-7) 

Next  we  try  to  find  the  expression  for  p(rfc|sfc,s/)jb,aj). 

At  time  k ,  the  probability  of  rfc  conditioned  on  sk,  sitk,  c,  and  C/  is 


p(l*fc|Sk,  &I,k>  C,  Cj) 

=  7r_L|R„|_1  exp  (rk  -  \fPacsk  -  \fPiCjSiikj  R“ 

(rk  -  y/Pacsk  -  S/Picrsj,)  | 

,  f  I  rke-y/PsCtSk-  y/PrCiyeSitk'2 

-  exp  <  —  2_j - 

l  1=1 


=  »-LiR«r 


at 


(6.8) 

(6.9) 


Assume  sk,siik  and  nk  are  independent  in  the  time  domain,  the  conditional  probability  of 
all  the  received  signal  rk  would  be 

Pfofclsk»  S/,fcj  ci  c/)  . 


—  n  P(r*l^*>  sI,h  ®j) 

i=k-{K- 1) 


=  7 r  LK\Ttn\  ^  exp  ^ 


E  E 

^  i=k-(K-l)  1=1 


n,t-VP3CtSi  -  s/PjcitSlj 


Denoting  the  exponent  in  (6.10)  as  A,  and  expanding  it, 
A  4  - 


VFsCiSi  \/PiCi^Si 

2 L,  2s  ^ 

i=k-(K-l)  1=1  1 

L  „  k  L 


-Ei  E  \n/-Ei[KP'ah-'/Pivi^- 

e=i  1  i=k-(K- 1)  t=i  1 

y/Pimcle]  [KPsa$+  (y/pfivtijt  ~  V^Ve)  ct  + 

{\fP$PiV* ci ,e  —  y/PsVtj  cjj , 


(6.10) 


(6.11) 
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where 


vi  =  e  r«*;  (6-i2) 

<=*-(*■- 1) 

=  E  ryS;,(  (6.13) 

i=fc— (if— 1) 

„  =  e  (6-14> 

i=k— (if — 1) 


Note  that  yt,yi,t  and  y  are  functions  of  s*  and/or  s 
Define 


Co  =  Ei  E  m* 

(6.15) 

/=i  *  i=fc-(if-i) 

Ci  =  ~  ~  \/^2//,<c/^] 

(6.16) 

x,  =  |x*|  e*^  =  \fP^Pjyc*It  —  y/PsUti 

(6.17) 

then  (6.11)  becomes 

L  l  1 

A  =  -C0  -  Y)  Ct  -  Y]  ~2  {KPsCc}  +  2\xt\  at  cos  (&  +  tpt)] . 

(6.18) 

t= i  *=i  ^ 

Substitute  it  into  (6.10),  then 

pfofelSfc.Sj^C.Cj) 

(6.19) 

(  L  N 

=  Ciexp|-]T]C*j 

(6.20) 

exp  |  ^2  +  2  |(E*|  Q!£  COS  (0£  +  i[>t)\  | , 

(6.21) 

where 

C\  =  tt~lk\Rz\~k  exp  {-C0}  • 

(6.22) 

As  mentioned  above,  the  channel  c  is  assumed  unknown.  We  eliminate  c  in  p(rfc 
by  integrating  over  it, 

Isfc)  S/,fc>  CI  ci) 

pfefclsfc)  ci) 

=  J P(rfc|sfc,  s /)fc,  c,  C/)pc  (c)  dc 

=  J  j  p(rk |sfc,  s/>fc,  c,  Cj)pa  (a)  p*  (</>)  dadcf), 

(6.23) 
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where  pc  (c) ,  pa  (a)  and  p#  (<£)  are  the  probability  density  functions  of  c,  a  and  (f>  respectively. 
Since  the  channels  are  independent  to  each  other, 


L 

Pa  (a)  P*  (<t>)  =  (6-24) 

e=i 

where  pae  (af)  and  p^  (<pe)  are  the  probability  density  functions  of  at  and  <f>t  respectively. 
For  Rayleigh  fading  channels, 

Pa((<xt)  =  ^exp(_§)  0<^<o°  (6-25) 

p<h(<t>l)  =  0<<fo<27T,  (6.26) 

where  Qt  =  E  [af]  is  the  mean  square  of  the  amplitude  of  the  ^th  channel. 

Substitute  (6.24),  (6.25)  and  (6.26)  into  (6.23),  after  some  straightforward  manipulations, 


p(r*lsfc>s/,*>c/) 


=  c>n 


\xtf 


-ct 


(6.27) 


jj  (KPsQe  +  a2)  j 
Expanding  the  exponent  in  (6.27)  and  collecting  the  terms  that  contain  aiti  and  <£/,*, 

VLt\xi\2 


a\  (KPsQe  +  crj)  °l 
=  Dt,i - i  +  2  \Dti3\ cxj cos  (<^/,<  ipt)  > 


(6.28) 


where 


n  -  g^M!! 

e>l  o\  ( KPSlt  +  al) 

_  QtPsPjjyl 

l’2  a\  ( KPSh  +  a}) 

Dt,  3  =  \Dt>3\e^  __ 

(KPSQ(  +  al)  s/Pi 
aj{KPsQt  +  aj)  m 

Substitute  (6.28)  into  (6.27), 


$hPWPj_ 

a2((KPsne  +  aj)yye' 


(6.29) 

(6.30) 

(6.31) 


P(l*lsfc>s/,k>c/) 

C^KPsne  +  a2 

exp  i  D(ti  —  +  DtfiOj t  +  2  |D^3|  ajjt cos  (<£/,*  —  <pt) 

I  o$ 


(6.32) 
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In  (6.32),  the  probability  of  r*  is  dependent  on  (£  =  1,  •  •  •  ,  L)  (which  is  included  in  Cj). 
We  can  cancel  this  dependence  by  integration, 

p{LkK,Si,k,ai)  =  j pfob |sfc,  S/,fc,  Cj)p^  (<t>i)  (6.33) 

where  p<p,  (<£/)  is  the  probability  density  function  of  <j>i.  For  independent  Rayleigh  fading 
channels, 

P<h  ( fa )  =  (^:)  0  <  fa, i> -  •  *  >  fa,L  <  27r  (6.34) 

Substitute  (6.32)  and  (6.34)  into  (6.33)  and  carry  out  the  integration,  we  have 


f>fofek,S/)fc,0!/) 

=  C, 


.{n 


Vlsdd  (sfc,  Si,k)  >  (6.35) 


where 

L 

rjMSDD  (s k,  sI  k)  =  exp  {£>£,1  +  h  (2  l-D^I  «/,* )  >  (6.36) 

£=1 

and  I0(x)  is  the  zeroth  order  modified  Bessel  function  of  the  first  kind.  Note  that  i)msdd  (s*,  sJifc) 
is  a  function  of  s*  and  sjtk. 

In  (6.35),  only  77 msdd  (sfc,Sj^)  dependent  on  (s/c,S/  ^) .  Maximizing  |sj.,s  j}k,ckj) 
with  respect  to  (sfc,  s/  fc)  is  equivalent  to  maximizing r? msdd  (sfc,  sI>k)  .  Therefore  tjmsdd  (s*,  sI  fc) 
can  be  used  as  decision  statistic  for  the  MSDD  detector.  The  corresponding  decision  rule  is 

(sk, s/,fc)  =  arg  max  Vmsdd  (s *,  sI<k)  .  (6.37) 

Sfc)8jifc 


The  MSDD  detector  searches  through  all  possible  (sfc,s/fc)  and  chooses  the  pair  that  has 

the  largest  tjmsdd  (s*,  s/>fc)  as  the  detected  output. 

To  complete  this  section,  we  briefly  review  maximum  likelihood  (ML)  detector.  It  is 
used  in  the  simulation  results  section  for  comparison  with  MSDD.  ML  detector  is  a  kind  of 
coherent  detection  technique  that  requires  the  channel  gain  of  both  the  desired  signal  and 
interference.  It  makes  symbol-by-symbol  detection  instead  of  sequence  detection  for  MSDD. 
The  ML  decision  rule  is  given  by 

(sfc.s/,*)  =  arg  max  p(rfc|s*,sj,*,c,c,).  (6.38) 

Sk ,si,k 

From  (6.38)  and  p(rfc|sfc,  s/,*,  c,  Cj)  shown  in  (6.9),  we  can  get  the  equivalent  ML  decision 
rule  as 

(sk,  si,k )  =  arg  max  t}Ml  ( sk ,  Si,k ) ,  (6-39) 

Sk,Sl,k 

where  the  decision  statistic 

1 2 

/  X  v1'  I rkt-y/P^CtSk  -  \/PicittSitk\  ,  . 

Vml  ( Sk ,  si,k)  =  2_s - ^2  • 

£=1  e 
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6.4  Simulation  Results 


The  communication  systems  simulated  had  4  diversity  branches.  DPSK  modulation  was 
employed.  The  channels  of  both  the  desired  signal  and  the  interference  were  assumed  to  be 
uniform,  i.e.,  =  1,  E  [o^]  =  1  and  of  =  a1  for  i  =  1, 2,  •  •  •  ,  L. 

In  the  figures,  SNR  is  the  signal  to  noise  ratio  defined  as  SNR  =  Pa/&2,  while  SIR  is  the 
signal  to  interference  ratio  defined  as  SIR  =  Pa/ Pj.  Fig.  6.1  and  Fig.  6.2  were  generated  for 
SIR  =  10  dB  and  SIR  =  -10  dB,  respectively. 

The  curves  labeled  “MSDD(tf  =  2)”  and  “MSDD(RT  =  7)”  are  the  results  for  the  MSDD 
detector  developed  in  this  chapter.  It  is  observed  that  the  performance  improves  with  the 
increase  in  K,  which  is  the  number  of  symbols  in  the  observation  interval.  For  example,  in 
Fig.  6.2,  at  BEP  =  2  x  10-3,  the  required  SNR  for  K  =  2  is  about  8.5  dB;  for  K  =  7  it  is  5.5 
dB.  That  means  increasing  the  observation  interval  from  K  =  2  to  K  =  7  symbols  results 
in  a  3  dB  SNR  improvement. 

The  curves  labeled  “OC”  are  the  results  for  OC.  The  curves  labeled  “MSDD(known  cov, 
K  =  13)”  are  for  the  MSDD  detector  discussed  in  Chapter  5,  which  was  developed  for 
known  covariance  matrix  of  the  interference  plus  noise.  MSDD  ( K  =  7)  has  about  the  same 
computation  complexity  as  MSDD  (known  cov,  K  =  13). 

In  Fig.  6.1,  the  BEP  of  MSDD  ( K  =  7)  is  larger  than  that  of  MSDD  (known  cov,  K  =  13) 
and  OC.  In,  Fig.  6.2,  the  BEP  of  MSDD  ( K  =  7)  is  less  than  that  of  MSDD  (known  cov, 
K  =  13)  and  OC.  We  conclude  that  at  high  interference  level,  MSDD  (K  =  7)  has  better 
performance  than  MSDD  (known  cov,  K  —  13)  and  OC.  That  can  be  explained  as  following. 
MSDD  ( K  =  7)  detects  the  interference  signal  s as  well  as  the  desired  signal  s*,  but  the 
MSDD  (known  cov,  K  =  13)  and  the  OC  detector  only  detect  the  desired  signal.  MSDD 
developed  in  this  chapter  is  a  kind  of  multiuser  detection  that  gets  better  performance  with 
the  increase  in  interference  power.  These  results  are  reflected  in  Fig.  6.3,  which  shows  the 
difference  of  the  required  SNR  of  MSDD  and  OC  at  BEP=  10  3. 

In  both  Fig.  6.1  and  Fig.  6.2,  the  performance  of  MSDD  ( K  =  7)  is  not  as  good  as 
that  of  maximum  likelihood  detector  (curves  labeled  “ML”).  But  the  difference  is  small  for 
a  high  interference  level.  In  Fig.  6.2,  at  BEP  =  2  x  10~4,  the  difference  of  the  required 
SNR  is  about  1.6  dB.  These  results  are  reflected  in  Fig.  6.4,  which  shows  the  difference  of 
the  required  SNR  of  MSDD  and  ML  at  BEP=  10-3.  We  can  expect  that  this  difference  will 
decrease  with  the  increase  of  K. 
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Figure  6.1:  BEP  versus  SNR  for  MSDD,  OC  and  ML.  For  L  =  4  branches,  SIR  =  10  dB. 
’MSDD  (K  =  2)’  and  ’MSDD  (K  =  7)’  are  BEP  for  the  MSDD  developed  in  this  chapter, 
while  ’MSDD  (known  cov,  K  =  13)’  is  for  MSDD  derived  in  Chapter  5 


Figure  6.2:  BEP  versus  SNR  for  MSDD,  OC  and  ML.  For  L  =  4  branches,  SIR  =  -10  dB. 
’MSDD  (K  =  2)’  and  ’MSDD  {K  =  7)’  are  BEP  for  the  MSDD  developed  in  this  chapter, 
while  ’MSDD  (known  cov,  K  =  13)’  is  for  MSDD  derived  in  Chapter  5 
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Chapter  7 

Summary  and  future  work 


7.1  Summary 

In  this  report,  we  present  the  following  work  we  have  done: 

•  Obtained  closed-form  expressions  of  the  exact  BEP  of  OC  for  BPSK,  multiple  inter¬ 
feres,  with  the  number  of  interferers  less  than  the  number  of  reception  branches. 

•  Formulated  simpler  asymptotic  expressions  of  BEP  of  OC  for  M-PSK,  one  interferer. 

•  Developed  the  decision  statistic  of  MSDD  for  communication  system  with  one  inter¬ 
ferer.  The  performance  of  this  detection  scheme  was  analyzed.  Through  analysis 
results  and  simulation  results,  we  proved  that  with  an  increasing  observation  interval, 
the  performance  of  MSDD  approached  that  of  OC  with  differential  encoding. 

•  Evaluated  the  performance  of  MSDD  for  the  case  when  the  channel  information  of  the 
interference  was  known. 


7.2  Future  Work 

Some  future  work  includes: 

•  Derive  the  closed-form  expression  of  the  exact  BEP  of  OC  for  BPSK,  multiple  inter¬ 
ferers,  with  the  number  of  interferers  is  equal  to  or  greater  than  that  of  the  reception 
branches. 

•  Try  to  derive  closed-form  expressions  of  the  exact  BEP  of  OC  for  M-PSK,  multiple 
interferers. 

•  Derive  simpler  approximate  expressions  of  BEP  for  OC. 

•  Apply  the  above  MSDD  scheme  (for  known  covariance  matrix  of  interference-plus-noise 
matrix)  to  multiple  interferers.  Analyze  and  compare  its  performance  with  that  of  OC. 
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Appendix  A 

Derivation  of  the  Characteristic 
Function  for  OC 


In  this  appendix,  we  derive  the  expression  (3.14)  for  the  characteristic  function  $d(jw)  of 
the  test  statistic  D. 

Define 

dm  =  A m9mX m  +  Am  Qmxmi  (A-l) 


then  from  (3.13)  we  have 


D  =  '£dm.  (A-2) 

m=l 

From  the  signal  model  in  Section  3.2,  the  definition  of  the  whitened  interference-plus- 
noise  vector  x  and  modified  channel  vector  g  in  (3.11)  and  (3.12),  after  some  algebra,  the 
covariance  matrix  of  x  and  g  can  be  evaluated  as 

R**  =  E  [xxff]  =  PslL  +  A  *  (A.3) 

Rgg  =  I L  (A-4) 

Rxg  =  R9X  =  y/P,lL.  (A.5) 

To  use  the  results  in  [7,  Appendix  B],  we  identify  the  following  quantities  using  the 
notation  in  the  reference:  Xm  =  xm,  Ym  =  gm.  Then  from  (A.3)  to  (A.5),  we  have  in  the 
notation  of  the  reference 


o 

II 

II 

1* 

(A.6) 

^xx,m  —  2  ^m) 

(A.7) 

1 

Myy.m  —  ^ 

(A.8) 

M xy,m  =  f^yx,m  =  ^ 

(A.9) 
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Substitute  the  above  equations  into  Eq.  (B-6)  and  Eq.  (B-5)  in  [7,  Appendix  B],  together 
with  Am  =  Bm  =  0  and  Cm  =  A"1.  After  some  straightforward  manipulations,  we  get  the 
characteristic  function  of  dm  as 


4>dm  (ju)  = 


_ 1 _ 


(A.10) 


where 

pi  ,171  =  V^T+  A  m  —  y/Ps  (A.  11) 

=  \fPa  +  Am  +  \/~Pa-  (A.12) 


It  follows  that  the  characteristic  function  of  D  is 

L 

=  n  **•  W") 

i 


(A.13) 

(A.14) 


m=l 

L 


=  n 


A  ('--Kir)  (1+^“ir) 


n 


(A.15) 


[i  -  jw  l1  (a?  ~  + 

For  ease  to  manipulate  later,  we  define  nm  as  (  notice  that  Am  =  Nq  for  m  —  Nj  + 1,  Nj + 


Mm  —  S 


+  a[K  + 1 

Am  ^  V  ^  A” 

_^5__ 


Ps_  _  /ZS 

-Af/  V  Am-. 
Moo 
A*oi 


+  ^ — — 
ATj  Am-Nj 


m  =  1,  •  •  •  ,iVj 

m  =  JVj  +  1,  •  •  •  ,  2  Nj 

m  =  2A/  +  1,  •  •  •  ,27V/  +  (L  -  A,) 
m  =  2NI  +  (L-NI)  +  l,---  ,2L 


(A.16) 


where 


Moo  = 


Moi 


iV0  V  ^ 


An 


V  A?  +  An' 


(A.  17) 
(A.  18) 


Then  $o(jw)  in  (A.15)  could  be  expressed  as 


2L  1 

= n  r^T 


m=l 


(A.19) 
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Appendix  B 

Evaluation  of  the  Residues  for  OC 


In  this  appendix,  we  evaluate  the  residue  in  (3.22). 

According  to  the  complex  variables  theory  [38] ,  if  a  function  /  (w)  has  a  pole  Wo  of  order 
N,  the  residue  of  /  (u;)  at  uiq  is 

Res  u  m  ;  u*]  =  (N  1 1}!  [(«  -  <*f  f  (“)]  L  •  (B1) 


*n  (3-21)  can  be  expressed  as 


(B.2) 
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B.l  Evaluate  the  residues  at  poles  u  =  — (m  =  Nj  4- 1,  •  •  •  ,  2 Nj 


For  Nj  +  1  <  m  <  2iV/,  the  poles  a;  =  —  j—  are  of  order  1. 


Res 


.  1 
J  J 

U  /V 


L_/ 

-A'll 

$!>(?£*>) 

r  v 

Mm/. 

a; 

0)=  — i-i- 

17  Mm 

i 

/2V/  A 

Wl*S 

* 

1  m_1  1 

1  TT  1 

H)2L 

(MooMoi)L  Ni  “  „AJf  (u  +  j+.'j 

2Ni 

TT  - 

1 

1  1 

11  / 

n=m+l  1 

'w+*i) 

(  .  i  \*-W/  /  , 

(w  +  '£)  ("  +  j£j 

u}=—j~L~ 
J  Mm 


(-i) 

2Nj 

n 


L_  fr\  1_\  1  1  IT  1 

j)2L  \n=l  f^J  (M 00Moi)L  n=l 

1 


n-m+1  (  +  j-^j  +  j^) 

After  some  simplification,  we  have 


L-V/  ' 


H«fhW  ji 

.  W  MmJ 


(^00  -  (rt.i  -  ' 


(B.3) 


(B.4) 
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B.2  Evaluate  the  residue  at  pole  u  =  —j 


Since  the  pole  co  =  —jj-^  is  of  order  ( L  —  V/), 


co  Moi 


dL~N >-x  j  r  (  .  i 

7,.,L-Nr-1  1  W  l  ^  i in- 


(. L-Nj  -  1)! 


-i —  fff  — ) - - 

■jfL  /**»  /  (^00/^01 


(/iooA«oi)L“iV/  (L-Nj-  1)! 


dL-Nl-l  ]_ 


du)L~Ni- 1  I  a; 


Define  /xq  =  oo,  then 


n=l(^  +  ii)J  (w+ii) 


where 


For  n  =  0, 


n— 1  J  2Nj  J 

n  (w  +  i; i)  J0L  (w  +  -?i) 

n— 1  2W/ 

_  /  i  \^Tr  tt  EiEn  TT 


2iVj  2Nj 

A0 = (-i)N/  tt  ,  ^ = (-i)^  ^since  = °°)  • 

££(/«> -/*«)  it 


For  1  <  n  <  2iV/, 


n-l  2N/ 

,  /  ■,\NrTl  ViVn  TT  W/i" 

*  =  (_1)  0(^y  n,(^ 

2  AT/  n— 1  t  2 


z/vj  n— i.  1  1 
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Using  (B.6)  and  the  derivation  expressions 


we 


dn 


duj 


-[/  MsM] 


SC) 


n\  dn  ef  (w)  dlg  (w) 

£  J  dun~e  doje 


<P 


1 


dun  (w  +  a)m 

get  the  derivation  in  (B.5)  as 


=  (-1) 


n  (ro  +  n-1)!  1 

(m  —  1)!  (u  +  a)m+n’ 


(B.10) 

(B.11) 


- > 

l 

£ 

1 

1. 

"  2Nj 

TT- 

1 

1“ 

11 

n=l 

L-Nj- 1  , 

-  E 
*=0  v 

L—Nj—1  , 

=  E 

t=o  v 


,  L-Ni 


Nj  -l 

\  dL-N’~ l~l 

~2N/ 

V- 

A 

de  1 

£ 

J  du1-”'-1-* 

n=0 

(W+J'i). 

duj1  (  .  i  \i_Ar/ 

(w+w 

2Nj 


L-  Nr- 
£ 

An 


L-Ni-l-l  (1  +  £  ~  ~  1  ~  l  ~  ^)- 

(1  - 1)! 


n=0  (w  +  jl) 

1 

After  some  simple  algebra, 


l+L-ATj-l-* 


(-1) 


(L-NT-  1)! 


(B.12) 


dL-N,-l 

J 

f. 

~2Nj 

TT 

1 

1 

dujL-Nj-1 

{“ 

11 

n=l 

("+*£). 

("+*») 

L-Nj 


=  (-1) 


L-Nr-l 


L-Nt-  1 


£=0 


(L  —  Ni  +  £—  1)! 


2Nj 

E 


A, 


[n=°  (w+i^) 


L—Nj—t 


(“+i  i)1-"'"'' 


(B.13) 
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(B.14) 

(B.15) 


*c(jw)  1 

>  3 

u  Hoi 


1 _ 1  xL-JV/-1  L  'ST'  1 

tom)L-Nt(L-Ni-ir  }  U 


(L  —  Ni  +  £  —  1)! 
t\ 


2  JV/  ,  \L—Nj—l  fn- 1 

L—Nj—t  ,  V'  _ „2JV,-1  I  TT _ - 

%  +  2^  (  xL-Nj-t^n  l  11  ...  _ 

^  (Hn  ~  MOl)  \i=l  ^ 


5*-*  L5iw-^ 


(/Xqq//qi)L  N/+< 


(B.16) 
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Appendix  C 


Express  Fj_n ,-k  (Pm)  as  a  Summation 
of  (yjm  —  rf)  to  Integer  Power 

In  this  appendix,  we  derive  (3.50),  which  expresses  the  function  Fk  (ym)  in  (3.45)  as  a  sum¬ 
mation  of  ( y ^  —  rf)  to  integer  power. 

For  simplicity,  define 

Fk  ~  (id*  Pm)  (jl  Vm)  d"  (1  Pm)  (jl  d~  2/m)  •  (C.l) 

Then 

F0  =  -22/m  (C.2) 

Fi  =  2  (1  —  7/)  ym.  (C.3) 

By  substituting  in  F*_ 1  and  Ffc_2,  it  could  be  easily  proved  that 

Fk  =  2r jF^  +  (2 L  ~  V2)  Fk- 2.  (C.4) 

To  simplify  the  notations,  define  mathematical  symbols  P  =  2t],  Q  =  y^—if.  Then  (C.4) 
becomes 

Fk  =  PFk-i  4-  QFk-2-  (C.5) 

We  try  to  find  the  relation  between  Fk  and  F0,  First  we  express  Fk  as  a  function  of 
Fk-2  and  F*_3.  From  (C.5), 

Fk  =  PFk-i  + QFk-2 

=  P  (PFk-2  d-  QFk-3)  +  QFk-2 

=  (P*  +  Q)Fk-2  +  PQFk- 3  (C.6) 

Continue  on  in  this  way,  we  have 

Fk  =  (P3  +  2PQ)Fk-3+(P2Q  +  Q2)Fk-4  (C.7) 

Fk  =  (P4  +  2P2Q  +  P2Q  +  Q2)  Fk-i  +  (PZQ  +  2PQ2)  Ffe_5.  (C.8) 

Judging  from  (C.5)  to  (C.8),  we  guess  that  the  relation  between  Fk  and  Fk-i,  Fk-r- 1  (/ 
is  any  positive  integer)  is: 
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•  when  I  is  odd 
o, 


*=[g(7;y-’v 

where  0/  =  (7  —  1)  /2. 

•  when  7  is  even, 

r 


Fk-i  + 


2tgt 


Pfc-/  + 


O/ 

E 

Lt=o 


/  - 1  - 1 
t 


yI-2t-lQt+l 


Fk-i- 1.  (C.9) 


P*-/-i-  (C.10) 


Lt=o 

where  0/  =  7/2. 

(C.5)  to  (C.8)  have  shown  that  (C.9)  and  (C.10)  are  valid  for  7  =  1,2, 3, 4.  The  validity 
of  (C.9)  and  (C.10)  for  any  positive  integer  7  could  be  proved  easily  by  the  method  of 
mathematical  induction. 

From  (C.9)  and  (C.10),  we  can  get  the  relation  between  7*  and  Pl,  P0  as  follow: 

•  When  k  is  even 

Let  I  =  k  —  1  and  substitute  it  into  (C.9), 


f „  =  |e  (*  j  ~ ‘V-'-vl  F>  +  [e  (* ” 1 7 1 " ()p*-i-m-iq*+i 

,t=o  \  t  /  J  Lt=o  '  ' 


F0.  (C.ll) 


Define  T  =  k/2.  since  I  =  20 j  +  1, 


°' =  ~T~  =  =  ^~1  =  T~1'  • 

Substituting  (C.12)  into  (C.ll),  we  get  the  relation  between  Fk  and  Fi,  P0  as 


(C.12) 


Ft=  J2(k  1  t)pi-1-2tQ‘  K+  £)(*  11  ()p‘-i-»-iQ‘+i 

.(=0  \  t  J  t_0  V  / 


To.  (C.13) 


•  When  k  is  odd 
Similarly,  we  can  get 

T 


Fir  = 


2£  ✓")£+! 


Ql 


F0  (C.14) 


where  T  =  (k  —  1)  /2. 

Finally,  substituting  P  =  2r?,  Q  =  -  ??2,  Po  =  -2ym,  and  Pi  =  2  (1  -  77)  ym  in  (C.13) 

and  (C.14),  after  some  manipulations,  we  have 

[k/2] 

Fk  (ym)  =  2 ym  ^  °*>*  ~  ^72)*  (c-15) 


t= 0 


where  [ik/2]  denotes  the  largest  integer  that  is  less  than  k/ 2,  and  ak,t  is  calculated  differently 
for  when  k  is  even  or  odd  as  follow: 
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•  When  k  is  even, 


(2v)*‘l  (1  -  V)  t  =  0 

[(‘T')  (I-*?)  -  2»(‘7ir*)]  1  <  <  <  [fc/2]  -  1 

-1  t  =  [*/2] 


(C.16) 


•  When  k  is  odd, 


(277)*"1  (1-7?)  t  =  0 

[fr)(l-r7)-  277(--)](277)fc-1-2‘  1  <  t  <  [*/2]  * 


(C.17) 


If  we  assume  (™)  =  0  for  m  <  n  or  n  <  0,  we  can  express  (C.16)  and  (C.17)  as  a  single 


expression: 


°‘.<=[C  1 0<i-’,)-2’'(W). 


\fe— 1— 2t 


for  0  <  t  <  [kj 2] . 
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Appendix  D 

Evaluation  of  the  Sum  of  Integrals 


In  this  appendix,  we  evaluate  the  sum  of  integrals  f  fm  (y)py{y)dy.  The  basic  steps  are: 

1.  Combine  the  sum  of  Nj  integrals  (each  integral  is  a  Nj— fold  integral)  into  one  Nj— fold 
integral. 

2.  Change  the  integration  limits  of  some  variables  from  finite  to  semi-infinite. 

3.  For  TV/  >  1,  separate  the  A/ -fold  integration  into  NT  independent  integrations. 


D.l  Combine  the  sum  of  TV/  integrals 


By  expressing  the  integrals  into  a  slightly  different  form,  we  can  convert  the  sum  of  Aj 
integrals  into  one  integral.  We  first  consider  the  integrals  /  fm  (y)  Py(y)dy  for  1  <m  <  Nr, 
then  for  m  =  1  and  m  =  A/. 

For  1  <m<  Aj,  /  fm  (y  )py(y)dy  is  an  A/- fold  integral.  We  first  carry  out  the  integra¬ 
tion  over  ym,  then  over  yNn  next  yNj~ i,  at  last  over  yx.  Since  oo  >  yx  ^  y2  ^  ^  Vn,  ^  V, 

f  fm  {y)py(y)dy  can  be  expressed  as 


/  fm(y)py(y)dy 

=  /  U  U  U  -{/  {L 


■  dym+ 1}  dym~i}  ■  •  •  dy2}  dyu 


(D.l) 


with  the  inner  integrations  being  carried  out  before  the  outer  integrations. 

From  (3.43)  and  (3.46), 

f  (  \  f  I  1  Vm  (1-*72)L  N‘  /tT  1  2/^_  \  /  TT  1  ~  Vn  \ 

fm  (y)Py(y)  -  2ym  {yl  _  ^  j  11  |nll+i  yl-yl] 


Kx 


Ni 


JJexp  [-/3  (y2  -  y2)]  (y2  -  y2) 
2\2 


2\  L—Nj 


.i=l 


n  &  -  yf> 

Li<*<j<W/-i 


ym-'-yNr 


(D.2) 
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Convert  the  variables  from  y  to  t  as: 

yi  -+  <1, 2/2  ->■  h  •  •  •  , 2/m— 1  ->■  tm-u  Vrn+ 1  tm,  ■  ■  ■  ,yN,  ^  tNj- 1,  ym  -4  .  (D.3) 

Then  (D.l)  becomes 

J  fm(y)py(y)dy 


ru'-iru 

•  *  *  rfifn}  ^m-l}  *  *  *  dt? } 


tm—2  (  ftm-1  (  rtNj-2  f  ftm-1  ~] 

[J  |  Jt  fm  (y)  Py(y)dtNj  |  dtNl- 1 1 


and  (D.2)  becomes 

1  +  (\  rfi\L~Nl  ( Nl 

r  f  \  /  \  1  -  tNj  (1  -  V  )  J  T 

/m(y)Py(y)  =  o7 - 7^ - ITl^F 

2^/  (4,  -  V2) 

[Nj 

Kx 


V/-1  1  _  /2 

n=l  %  M 


nexpb^-V)]  (t*  —  V2)L  Nl 

2  +2''2^  tit2  ■  ■  '  tffj, 


.*=1 


n  «-« 

Ll<i<i<AT/-l 

which  could  be  changed  to  (by  re-arranging  the  terms) 

/m(y)  Py(y) 

=  iKr  (1  -  ri‘)L-m  (n  (1  -  «  exp  l-fi  «  -  p2)]  (t2  -  rf)1-"’  <» 


t  n=l 
2  *2\2 


n  « - 

Li<i<j<V/-i 

n‘(4-« 

.  n=l 

Substitute  (D.6)  in  (D.4), 


(1  -  tNt)  exp  [-p  (4,  -  t?2)] 


J  fm(y)Py(y)dy 

=  iK^r-rf'-r 

&  Jrj  Jr]  Jr\ 

(n^l-^expf-^t2-,2)]  (tl-y*)L-Nlt, 

l  n=l 

Utm—1 

(1  -  tNl)  exp  [-p  (4,  -  V2)] 

dttfj—i  •  *  *  dt^dti. 


•i 


(D.4) 


(D.5) 


(D.6) 


n  w  -  <?)! 

U<i<i<V/-i 


Wj-1 

IT  (*N,  -  ^n) 

n=l 


) 


(D.7) 
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The  above  equation  is  for  1  <  m  <  Ni.  Similarly,  for  iVj  =  1  and  Nj,  we  have 

J  h  (y)py(y)dy 

"  «/T7  «/  7}  Jr} 


f  Ni  —  l 


n  c1  -  #)  exp  [-^  -  ^2)]  -  f?2)*1 1 


n=l 
roo 


*„} 


n  m-# 

X<i<j<NT~  1 


(1  -  fjv>)  exp  [-0  (4,  -  7?2)] 


AT/-! 


IK4.-4) 


n=l 


dt 


Ni 


dtpjj- 1  •  •  ■  dt2dt\ 


(D.8) 


and 


J  fNi(y)py(y)dy 

^  Jr]  Jr]  Jr] 

{Nj- 1 

n(l-4)expM«--!2)]«--!2)t_K'<. 

n=l 

{/  "  ^  ~  <7^exp  (4,  _  ^2)] 

•  dt2dti. 


AT/-1 

II  (4/  “  4n) 

n— 1 


4f  n  (4?-*?)2 

‘ ;  i 

'JV/  r 


(D.9) 


Notice  that  the  integrands  in  (D.7),  (D.8)  and  (D.9)  are  the  same.  The  integration  limits 
for  ttfj-i,  •••  ,h,ti  are  the  same  as  well.  Therefore  we  can  write  the  sum  of  the  JV>  integrals 
as 


X)  j  fm(y)py{y)dy 

m— 1  u 

"  «/r/  Jr?  Jr] 

In  (1  - Q exp l-p (t2 - 72)]  ft - n‘)L-Nl t, 


n=l 


<x 

'Ni- 1 

n  (4/  - t] 

Ln=l 


n}[  n 

J  Ll  <i<j<N/—l 

Nl~l  rtm- 1  ftNj- 1\ 

+  J2  I  +  f  )  (!  -  *n,)  exp  [-0  (t%  -  t?2)] 

m= 2  ^ ^  / 

n)|  | 


dtpij—i  *  •  *  dt%  dt\. 


(D.10) 
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Since 


X)  [  fm{y)py(y)dy 

m=lJ 

-  iMisr-rr-r 

-»  Jr}  Jrj  Jf] 

l  n= 1 


■) 


(1  -  tpfj )  exp  [-0  (t2Nl  -  r?2)] 


Nj-l 


n  -  *») 


n=l 


n  «-«’ 

dtffj  \ 


dttfj— i  *  •  *  dt^dti. 


(D  11) 


which  consists  only  one  AT/— fold  integral. 

To  simplify  the  notations,  we  convert  variables  as:  (t2  —  if)  — >■  zn  for  n  =  1, 2,  •  •  •  ,  N[—l, 
and  define 


20  =  1  -  »T, 


(D.12) 


then  (D.ll)  becomes 

Ni 


£  [  fm(y)py(y)dy 

-i  /*0O  fZl  fZNj-2 

-  &**!.  L  -l 

{ n  (z° ~ Zn ) exp zn~Ni j  n  <*- 

{r 


(1  -  *JV,)  exp  [-/?  (t2Nl  -  V2)] 
dzpjj— i  •  *  *  dz,2idz\ . 


n  r?2  2n) 


JV/ 


} 


(D.13) 


D.2  Change  the  integration  limits  of  some  variables 

Consider  the  integrations  in  (D.13).  The  integration  area  for  Z/y, -i,  ?  22,  2i  is  oo  >  z\  ^ 

22  ^  23  ^  ^  znj—i  ^  0.  There  are  other  similar  areas  such  as: 

OO  >  2l  ^  23  ^  22  ^  ^  ZNj—1  ^  0, 

OO  >  22  ^  2i  ^  Z3  ^  ^  ZNl- 1  ^  0, 
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etc.  There  are  (iV>  —  1)!  such  areas  in  the  region:  oo  >  zi,  z2,  zs,  ■  •  •  ,2^,-1  ^  0.  Since  the 
integrand  is  symmetric  to  z\,  z2,  •  •  •  ,  z^j-i,  the  integration  over  these  (iVj  —  1)!  areas  must 
be  equal.  Therefore  the  integration  over  one  such  area  is  1/  (iVj  —  1)!  of  the  integration  over 
the  region:  oo  >  zt,  z2,  z3,  •  ■  ■  ,Zf}j-\  ^  0.  Therefore 


Nl  r 

53  /  fm(y)py(y)dy 

m=l  J 


1 


:Ki$ 


L—Nj 


poo  poo  poo 

Jo  Jo  Jo 


(Nr  -  1)!  2Nj 

{ n  _  *•»)  exp  zn~Ni  1  n 

l  n=l  J  U<i<i<AT/-l 

Uoc  rw/-i 

(1  -  tNj)  exp  [-/?  ( t2Nj  -  rf)] 


(*  ~ 


n  (^/  ~r,2~ Zn) 

L  n=l 


dt 


N, 


dzffr- 1  •  •  •  dz2dzi. 

The  above  integration  is  treated  differently  for  Nj  =  1  and  Ni  >  1. 


CD-14) 


D.3  n,  =  1 

Prom  (3.48)  we  have 

(D.14)  simplifies  to 


Kx  = 


2  ql 


(L~  1) 


53  /  /m(y)py(y)<*y 

m- 1 17 

=  r  -  *>> exp  h*  w  -  ’,2)i  *>• 


Since  77  =  ^  +  1, 


^  =  1  -  r/2  = 


Using  the  above  equation  and  the  result  in  Appendix  E,  we  have 

Nj  1  /  i  \  L—  1 


tjf-  MMrMr  -  (~)  "  (B»  -  5?) 

m— 1  ^ 


where 


poo 

B0  =  /  (1  -  ii)  exp  [-0  (t?  -  t?2)] 

Jr) 

=  P  (/V)  Q  (y/Wv)  • 


(D.15) 

(D.16) 

(D.17) 

(D.18) 

(D.19) 
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D.4  Nj  >  1 

For  this  case,  we  first  try  to  separate  the  multiple  integration  over  Zut- 1,  •  •  •  ,  Z2,  z\.  Accord¬ 
ing  to  [39,  Chapter  3],  for  any  bounded  function  /  (•) , 


noo  roc  AT/-1 

•••/  n/(^) 

Jo  n=1 


n  -  *?■) 

U<»<j<W/-i 

1 


d^/-i  •  •  •  dzidzi 


,Ni- 1 


noo  roo  * 

n/w 

•/  0  n=1 

dzNl-i  •  •  •  dz^dzi 


Z2 


zl 


Nt-2  Ni- 1 

Z1  z2 


Nr— 2 
zNj- 1 
JV>-1 
zNi- 1 


y2AT/-4 

ZAT/-1 


(D.20) 


where  |-|  is  the  determinant  of  a  matrix.  The  element  on  the  ith  row,  j-th  column  of  the 
matrix  in  (D.21)  is  z^+j~2.  Notice  that  all  the  elements  on  the  j-th  column  of  the  matrix 
depend  only  on  variable  zj. 

Therefore  (D.14)  is  equivalent  to 


X)  [  fm(y)Py(y)dy 

m- 

1  poo  poo  pc 

L  -L 


{Nj- 1 

n 

n=l 


(20  -  zn)  exp  {-fiz, 


1 

Zl 


zl  z‘. 


(1  -  tNl)  exp  [-0  (fNt  -  V2)] 


‘Nj- 1 


L  n=l 


dzNj~  1  *  *  *  dz2dz\. 


ously, 


Z2 

Jit-  2 
ZNj~l 
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Z2 

JV/-1 

NT-l 

2 

*  *  *  A 

•'iVj-l 

\ 

-v2 

zn ) 

dtNj  ? 

) 

(D.21) 

n  ni] 

l 1 

- ri 2- 

2n) .  Obvi- 

~rn\ 

Z1 

(D.22) 

n=l 


P= 0 


mi  =0,1 

mi=l—  p 


n(*o-*,)=E(-i>2^  e  < 

p=0 


mi  m2 


n=l 


mi, m2  =0,1 
mi+m2=2- p 


(D.23) 
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In  (D.23),  the  second  summation  is  taken  over  all  the  sets  of  (mi,  m2)  which  satisfy  the 
two  restrictions:  (1)  each  m,  (i  =  1, 2)  must  be  0  or  1;  (2)  The  summation  of  all  m^s  must 
be  equal  to  2  —  p. 

Similarly,  we  have 


N,-l  Ni- 1 

f[(Z°-Zn)  =  5Z  (-1)^"1"Pz( 

n= 1  p=0 


mi  m2  m^-i 

Zl  Z2  ***^/-l 


(D.24) 


Nt-  1 


Nj-1 


n— 1 


q=0 


n\  ri2 
z2 


nNl-t 

zNi- 1 


ni+n2H - hnjv/-i=W/-l-<7 

".€{0,1} 


(D.25) 


Substitute  (D.24)  and  (D.25)  in  (D.21), 


X)  [  fm{y)py{y)dy 

m= 1  ^ 

AT/-1 

E 

p=0 


ATj-1 


-i-<7 


E 


fl=0 

E 


mi+m2H - (-mjVj-i=Af/— 1—  P  ni+n2+—+nNI-i=Ni~l-q 

mj6{0,l}  "<€{0,1} 

f*  oo  poo  poo 

ymi+ni  m2+ri2  .  #  .  mNj-i+nNj-i 


noo  roo 

■  z? 

Jo 


'AT,- 1 


n  eXP  (~PZn)  zn 


L- AT/ 


71=1 


1 

Zl 


&AT/- 1 


^2 


Nj—2  Nj-l 
zl  z2 


yNj-2 

ZNj- 1 
JV/-1 
1 


2AT/— 4 
^AT/-! 


{/  ^  ~  ^  ^  eXP  M  (&/  “  »J*)]  } 


'V 

dz^-i  •  •  •  dz^dzi 


(D.26) 


After  combining  2^ln+nnexp(— ^3zn)  (for  n  =  1, 2,  •  •  •  ,  Nj  —  1)  into  the  n-th  column 
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of  matrix,  (D.26)  becomes 


p. 

/  /m(y)py(y)dy 

m=l'' 

AT/-1  Ni- 1 

=  -L K e  e  (-1)'"-1 

^  -  g=0 


E  E 

mi+m2H - (-m/ifri=Wj-l-p  ni+ri2H — t-»iw7-i  — Ni—1  Q 

mie{0,l}  «{€{(), 1} 

Ui,i  Uifi  Ui'Nj-i 

U2,l  U2,2  •  ‘  ‘  U2,N,-\ 

UNj-  1,1  UNJ- 1,2  ••• 

(1  -  tNl)  (t2Nl  -  r?2)9  exp  [-P  (t2Nl  -  »72)]  dtjv;  j 
dztfj- 1  •  •  •  dz2dz\. 


where 

««  =  exp  (-0Z,)  tf"**™*™-*. 

Since  the  elements  on  the  j-th  column  of  the  matrix  only  depend  on  variable  Zj 
put  the  integration  over  Zj  into  the  matrix  and  obtain 


where 


Nl  r 

J2  /  fm(y)py(y)dy 

m=l  J 

ATj-1 

E 

p= 0 


Ni- 1 


9=0 


detV 


E  E 

mi+m2+,”+wijVj _i=iV/— 1— p  ni+«2H — t-njv7_i— AT/  1  9 
mje{0,l}  n<€{0,l} 

|  jT  (1  -  tNj)  (t%  -P2)q  exp  [-0  {t%  - 1 f) ]  dtNl | 


detV  = 


VM 

Vxa  ■ 

Vi.Ar.-i 

Vi,i 

V2,2 

V2,Nf-l 

Vat7-1,1 

Vnj-1,2 

•  •  *  Vv/— l.AT/— 1 
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(D.27) 

(D.28) 
we  could 


(D.29) 

(D.30) 


and 


roo 

T  r  /  /  Q  \  TTlj  +7T-J  +L— ATj+t+j“ 2  1 

Vij  =  /  exp  (-£*,•)  ^ J  <% 

J  0 


_  (mj  +  rtj  +  L  -  JV>  + 1  +  j  -  2)! 

^mj+nj +L- AT/ +i+j-2+l 

Substituting  (D.31)  in  (D.30)  and  taking  out  the  terms  with  /?,  we  have 

_ 1 _ 

detV  _  ^Efii_,K+nJ+i-N/+i-2+i)+Ejl/r1< 

Wi,i  Wl.2  •••  WitN,- 1 
^2,1  ^2,2  •  '  '  W^ATj-l 

WjVr— 1,1  W^iVj-1,2  "•  Wjv,— 1.AT/-1 

Wij  =  (rrij  +nj  +  L-NI  +  i  +  j-  2)!. 


(D.31) 


where 

Since 


(D.32) 


(D.33) 


Nj-l  N j- 1 

(m^  +  rij-  4-  L  —  Ni  +  j  —  2  + 1)  +  i 

j=i  *=1 

JV/-1  Ni-l  Nj—l 

=  +  j3ni  +  (L-^/)(iV/-1)+  Y/j-(NI-l)+ 

i=i  i=i  i=i 

=  Ni  —  1— p  +  Ni  —  1  —  q  +  (L  —  Ni)  (Ni  —  1) 


Nt-  1 


t=l 


+w-lm_w_1)+ 

=  (l  + 1)  (jv>  - 1)  -p-«, 


(iV/  —  1)  Nj 


detV  =■ 


•detW. 


'  p(L+l)(Nj-l)-p-q 

Next  let’s  evaluate  the  final  integration  in  (D.29).  Separating  the  first  term,  then 

roo 

I  (1  -  tNl)  (4,  -  r?) 9  exp  [-P  (t%  -  rj2)]  dtNl 
Jr] 
r°° 

=  /  (4,  -  7?2)?exP  [-P  (4/  -  v2)]  dtNl 

Jr] 

r°° 

-  I  tNl  (4,  -  v2)9  exp  [~p  (4/  -  n2)]  dtN, 

Jr] 

f°° 

=  I  (4,  -  *?T exp  [-0  (4/  ~  V2)]  dtN, 

Jr] 

-\j  (^-v2)9  exp  [-0  (4,  “  ^2) ]  d  (4,  “  ^2) 


(D.34) 

(D.35) 


_  o  I_iL 

9  2  09+1 


(D.36) 
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where 


roc 

Bq=  I  -  r?) 9 exp  [-P  (4  -V2)}  dtNl , 

J  TJ 


which  can  be  calculated  recurrently  as  discussed  in  Appendix  E. 
Substituting  (3.48),  (D.35)  and  (D.36)  in  (D.29),  we  get 


AT;  . 

53  /  fm(y)Py(y)dy 

171=  1  J 


PLN,ZL-H, 


[ni'i  (i  -  0!]  [n£.  w  -  <)! 

Nl- 1  JV/-1  /  i  n\  \ 


E 


mi-fm2H — hmifj-pAf;- 1—  p  tii+juH — Nj  — 1— g 
m<e{0,l}  nie{0,l} 

_ I _ Hpt  W 

p(L+l)(Nj-l)-p-q  061  VV- 


Define 


(D.37) 


(D.38) 


H  (p,  q) 

1 


n£,  (i  - ;)!]  [n£,  (n,  -  i>! 


E  E 

mi, m2,—  ^^-1=0,1  ni,n2,—  0,1 

mi  +m2+— +mAf7_i=Af/— 1— pni+n2+— +njV/-i=ATr— 1— 9 


detW. 


(D.39) 


Substituting  in  z0  =  and  after  some  manipulations,  we  can  simply  the  expression  in 
(D.38)  as 


53  [  fm(y)py(y)dy 

m=l  “ 

=  H) 


L-Nj  Nj—\  Nj- 1 

E  E  (-l)’H(p,q) 

p= 0  q— 0 


(D.40) 


which  is  the  expression  we  want  to  get. 

If  we  assume  H(p,q )  =  f°r  A/  =  1,  then  by  setting  TV/  =  1,  (D.40)  is  equal  to 

(D.18).  That  means  (D.40)  could  be  used  for  any  A/. 
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Appendix  E 

Derivation  of  Series  Bq 


In  this  appendix,  we  derive  the  method  to  calculate  the  series  Bq  (q  =  0, 1, 2,  •  •  •),  which  is 
defined  as 

poo 

B,=  I  exp  (E.1) 

In  the  following,  we  first  derive  the  expressions  for  Bq  and  B\.  Then  we  will  show  that 
for  q  >  2,  Bq  could  be  evaluated  by  Bg-i  and  Bg-2. 

E.l  q  =  0 


B0  =  j  exp  (— /?  (t2  —  rf))  dt  =  exp  (firf)  j  exp  {-pt2)dt.  (E.2) 

Define  \z2  =  fit2,  then  x  =  ^=jZ,  z  =  and 

B0  =  exp  {Pi f)  J  exp  {-fit2)  dt  =  exp  {Pi f)  exp  ~^pdz 

—  — ^exp  (firf)  V2kQ 

=  ^exP  (^?2)  Q  (^)  ■  (E*3) 


E.2  q  =  1 


poo 

Bx  =  f  exp  {-P  (l t 2  -  rf))  (t2  -  T]2)  dt 
Jrj 

r  poo  roc 

=  exp  (fiif)  [j  exp  {-Pt2)  t2dt  -  if  J  exp  (-0t2)dt 


(E.4) 
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E.3  q  >  2 

In  this  case,  we  express  Bq  as  a  function  of  Bq-i  and  Bq- 2.  The  details  are  as  follows: 

poo 

Bq=  exp  (-/?  (; t 2  -  r]2))  (t2  -  t/2)^1  (t2  -  r f)  dt 
Jrj 

=  r  exp  (-fi  (t2  -  ,’))  (t2  -  r?)'-'  t2dt 
Jq 

poo 

- r, 2  exp(-0(?-n2))(C‘-n2y-1dt 

Jr} 

/oo  1 

f  -  rf)q-1  — tdexp  (-/3  (t2  -  v2))  -  rfBq _i 


if  ~V2Y  1  P  (~P  f  -  V‘ 


1  f°°  r 

+  2p  4  GXP  ^  ~  7?2))  ^  f  ~  ^Y~2  2f2  +  (<2  -  ??2)9"1]  dt  -  fB<l- 

=  ^  1°° exP  H?  (i2  -  ^))  [2  (g  -  1)  if  -  V2)9-2  f  -rj2  +  rj2)  +  (t2  -  n2)*-1} 


=  2)3/  exp  (“^  “  ?>2))  [(2?  “  x)  (<*  -  *?T  1  +  2  (q  -  1)  t?2  (t2  -  V2)9  2] 

=  2/3  _  ^  59-1  +  2  ^  ^  -  V2Bq-l 


(2g  —  1) 

2/3 


.^>2  d  ,  1 )  „2  o 

—  V  I  "fl-l  3 - ~p — V  Bq_ 2. 


Therefore  S9  can  be  obtained  from  £L_j  and  £L_2. 
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Appendix  F 
Integration  of  Iq 


(F.2) 


In  this  appendix,  we  prove  the  following  integration  of  zeroth  order  modified  Bessel  function: 

f°°  1  B 

J  te~at2 10(2y/pt)dt  =  —  e«  (F.l) 

Eq.  (6.614.3)  in  [23]  is 

Let  v  =  0, 

r°°  s 

l 

Using  equation  of  M_/1_ili(x)  in  [40,  P.  432  ], 

r°°  ,  ! — ■„  e^i  ,B.i  ig.  1  i. 


(F-3) 


Let  x  =  t2, 


It  follows 


/°°  e-axI0(2\/]te)dx  =  [°°  2te~at2 Io (2 \f~0t)dt  =  - 

Jo  Jo  a 

pOO  1  Q 

J  te~at2 10(2y/j3t)dt  =  — e« 


1  & 
e° 


(F.4) 


(P-5) 
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Appendix  G 

Derivation  of  the  Characteristic 
Function  for  MSDD 


In  this  appendix,  we  derive  the  expression  (5.32)  for  the  characteristic  function  of 

the  MSDD  test  statistic  D.  To  that  end,  using  (5.24)  and  (5.30),  we  express  the  test  statistic 
D  in  quadratic  form: 


D  =  ^2  A\  (|y*(s*)|2  -  MsJt)|2)  » 


(G.l) 


where 


Xt  (KPS  +  Xi) 


for  £  =  1,2,--  •  ,L.  Define 


it  =  Aj  (|y,(st)|2  -  |jtt(syf)  , 


(G.2) 


(G.3) 


D-Z*. 


(0.4) 


Also  define  vector  y  (s*)  =  \yi  (sjt) , 3/2  (s it) ,  •  •  •  ,Vl  (sfc)]T  1  then  from  (5.12)  and  (5.8),  we 
have 


(Sfc)  —  [  ^  X-k-iSk-i 


=  uf  em. 


(G.5) 
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Prom  the  signal  model  in  Section  II,  E  [y  (s*)]  =  E  [y  (s*)]  —  0.  After  some  algebra,  the 
covariance  matrix  of  y  (sfc)  can  be  evaluated  as 


E  [y(sfc)y  (s,fc)ff]  =  E 


(K- 1 


(K- 1 


H 


U, 


,  »=o 


vf=0 


K2PsIl  +  KAz. 


Similarly,  we  have  the  following  results: 

B[yK)y(4)"]  =  |«(s*,s'i)|2P.iI  +  *rA, 

E  y  (s„)  y  (s»"j  =  t>*(s*X)  (h'P,lL  +  A,) 
y(4)y(s*)B]  =  w(s*>s*)  (epsil  +  a2)  , 


E 


(G.6) 

(G.7) 

(G-8) 

(G.9) 


where  v(s*,s'fc)  =  sjfs*. 

To  use  the  results  in  [7,  Appendix  B],  we  identify  the  following  quantities  using  the 
notation  in  the  reference:  Xi  =  ye(sk),  Ye  =  ye(s'k).  Then  from  (G.6)  to  (G.9),  we  have  in 
the  notation  of  the  reference 


X,  = 

V(  =  0 

(G.10) 

\  (K2P,  +  K\,) 

(G.ll) 

i(|^Sl,4)|2P.  +  A-A£) 

(G.12) 

Mxy,/  — 

io-^.slHifP.  +  A,) 

(G-13) 

llyx,l  — 

i„(St,sl)(i-P,  +  A,). 

(G.14) 

Substitute  the  above  equations  into  Eq.  (B-6)  and  Eq.  (B-5)  in  [7,  Appendix  B],  together 
with  A  =  Aj,  B  =  -A\,  and  Ct  =  0.  After  some  straightforward  manipulations,  we  get  the 
characteristic  function  of  de  as 


4>dt  (ju)  = 


Q\,t  @2  ,e 


(UJ  +  jd\/)  (w  —  j02,() 


(G.15) 


where 


**  =  (G16) 
=  t[KpJxSm  +  ca] .  (0.17) 


and 


C  =  tf2-Ms*X)l2. 


(G.18) 
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It  follows  that  the  characteristic  function  of  D  is 


<Mjw) = n  <t>dt  o'w) = n  y 


*=1 


t=l  ^1 


■*£)  (1+jl%) 


(G.19) 


For  a  system  with  a  single  interference  source,  the  eigenvalues  of  the  interference  plus 
noise  covariance  matrix  are  A t  =  cr2  for  l  =  2, 3,  •  •  •  ,  L.  It  follows  that  6\ti  =  #1,2  and 
02  t  =  02,2  for  i  =  2, 3,  •  •  •  ,  L.  Hence  the  characteristic  function  can  be  expressed  as 


$r>(jw)  ^  (1  _  jn2u>)  (1  _  oj)l  1  (1  -  jn iU))L  1  ’ 


(G.20) 


where  Hi  =  l/0i,i,  M2  =  —1/02, 1,  M3  =  1/01,2,  and  /x4  =  — 1/02,2- 
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